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WANL-TME-1906A 
NASA CR- 121 225 


AN UPDATE OF INPUT INSTRUCTIONS TO TEMOD 
I. INTRODUCTION 

The theory and operation of a Fortran IV computer code, designated as TEMOD*, 
used to calculate tubular thermoelectric generator performance is described in WANL-TME- 
1906. The original version of TEMOD was developed under AEC Contract AT(29-2)-2638 
in 1969. This report which is written as Appendix D of WANL-TME-1906, describes additions 
to the mathematical model and an update of the input instructions to the code which have been 
developed under AEC Subcontract N854-0051, in the period 1969 - 1973. 

Although the basic mathematical model described in WANL-TME-1906 has remained 
unchanged, a substantial number of input/output options have been added to allow completion 
of module performance parametrics as required in support of the Compact Thermoelectric 
Converter System Technology Program conducted at Westinghouse Astronuclear Laboratory. 

The report, then, basically replaces Section IV of WANL-TME-1 906 entitled "Input to the 
TEMOD Code." 

Section V of this report contains a Fortran listing of the code. 


* 


C. M. Rose, "A Numerical Model for Tubular Thermoelectric Generator Performance 
Analysis", WANL-TME-1906, April 1969. 



II. BASIC DESCRIPTION OF CODE SETUP 


The TEMOD code consists of a main program, referred to as TEMOD, seven sub- 
routines and three function subprograms. Although a complete Fortran listing of the code 
is given in Section V of this report, a brief discussion of the code sections is given below: 

A. MAIN PROGRAM: TEMOD 

Contains read statements for input quantities. Lists input parameters and module 
dimensions unless list is suppressed by input control parameter NZ(5). Also directs flow 
of logic between each of the individual subroutines and subprograms. 

B. BLOCK DATA SUBPROGRAM 

Contains compiled tables of all applicable material properties. Selections of 
thermoelectric material combinations, clad and conductor ring materials can be made 
by specification of control parameters at input. 

C. PHONY SUBROUTINE 

Assigns thermoelectric, clad and conductor ring properties as specified by input 
control parameters. Also adjusts thermoelectric material properties by percentages specified 
during input. Outputs all material properties used in each calculation unless output is 
suppressed by input control parameter NZ(5). 

D. SUBROUTINE COUPLE 

Performs heat balance/radial temperature profile/electrical output calculations 
for each thermoelectric couple. The mathematical model used for these calculations is 
discussed in WANL-TME-1906. A cross sectional view of a unit couple is shown in Figure 1. 

E. SUBROUTINE OPTIM 

Performs temperature, dimensional and load resistance incrementations specified. 
Also contains output statements which are restricted to one line of parameters for each 
set of conditions. 
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FIGURE 1. CROSS SECTIONAL VIEW OF A "UNIT" COUPLE 


F. 


SUBROUTINE PUMP 


Performs dimensional incrementations specified for pump module parametric 
calculations in which electrical power and required load voltage have been set. 

G. SUBROUTINE LIFE 

Performs performance as a function of time calculations in which operating condi- 
tions and degradation rates are specified. 

H. SUBROUTINE RITE 

Combines the results of the individual couple calculations to determine module 
performance for any specified number of couples (see NZ(10) below) operating in thermal 
parallel and electrical series. Also contains write statements to list results of module and 
individual couple calculations. 

I. FUNCTION SUBPROGRAM SI 

Performs all interpolation or extrapolation calculations. Primarily used to evaluate 

material properties which are contained in the program as temperature table values. 

J. FUNCTION SUBPROGRAM DK2FK 

Performs all temperature unit conversions. Performance calculations are done 
in Kelvin units, but input/output may be specified in either Fahrenheit or Kelvin units 
(see NZ(8)). 

K. FUNCTION SUBPROGRAM WATE 

Calculates module weight based on input module dimensions. Weight calculations 
will include contribution of end closures if appropriate entry is made in Z(61). 
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III. INPUT TO THE TEMOD CODE 

A. GENERAL 

There are four categories of input data required to operate the TEMOD code,, 

The formats used to read these parameters have not been modified from the description 
given in WANL-TME-1906. The basic categories of input data are: (1) bulk material 
properties, (2) fixed point (integer) control constants; (3) floating point data, and (4) oper- 
ating temperatures. 

B. BULK MATERIALS PROPERTIES 

Bulk materials properties include thermal conductivities and densities of all 
materials in the generator, Seebeck coefficients of all thermoelectric materials, and 
electrical resistivities of all materials in the electrical circuit. These properties with 
the exception of densities, are introduced as tables with the property evaluated at 50°K 
increments from 300°K up to 1000°K (15 values). 

A listing of the material property subroutine, DATAIH is given in Appendix A 
along with the complete program listing. Properties for six types of thermoelectric materials 
(TEGS-3Nr TEGS-2N or GE-nl, TEGS-3P, TEGS-2P, ternary n-type, and ternary p-type) 
have been built into the code and calculations can be made using any combination of 
these materials by proper selection of the control constant, discussed below. In addition 
properties are entered for three types of cladding materials (stainless steel, inconel and 
tantalum) and three types of conductor ring materials (iron, tungsten, and molybdenum) 
are also entered and can be selected by an input control constant. 

The sixteenth entry in the thermal conductivity table for each material in the 
DATAIN subroutine listed in Appendix A is the density of the material (in pounds per cubic 
inch). These densities are used in weight calculations. 

Also note that Seebeck coefficients for all thermoelectrical materials are entered 
using absolute values. The negative Seebeck coefficients of n-type materials is handled 
internally by the program logic. 
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c. 


FIXED POINT CONTROL CONSTANTS 


All fixed point data is read into a list called NZ. The list is dimensioned 50 
although not all of the 50 locations are used. This list below gives the instructions that 
correspond to each location in the NZ array. The method of entering numbers into this 
array is unchanged from the procedures specified in WANL-TME-1 906. 

Table I lists fixed point control parameter definitions for each entry in the NZ 
array. As shown in the table, the first four entries in the NZ array refer to thermoelectric 
materials which must be specified for the inner and outer radial segments of both the n- 
and p-legs of the thermoelectric couples. A schematic of a "unit couple" is shown in 
Figure 1. This option allows performance calculations for modules incorporating radially 
segmented thermoelectric washers. By specifying the same thermoelectric material to 
the inner and outer segment of either leg, results will correspond to modules in which 
no radial segmenting has been incorporated. 

As discussed above, property tables for six types of thermoelectric materials have 
been built into the code. By specifying an input control constant ranging from one to six, 
the material property tables for any of the six types of thermoelectric materials can be 
used in either radial segment of either leg. 

All entries designated by an asterisk in Table 1 refer to parameters which are 
zeroed at the beginning of each calcuiational case to avoid potentially expensive (in 
terms of computer time) errors. If non-zero entries are desired in change cases for any 
of these parameters, the values must be reset in each change case. All other entries in 
the NZ array will maintain their previous values for all change case, unless changed by 
entering new values in subsequent cases. 
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TABLE 1 

TEMOD INPUT FIXED POINT CONTROL PARAMETER DEFINITIONS 


NZ- = J: 


NZ(2) = J: 


NZ(3) = J: 

NZ(4) = J: 


NZ(5) = I RITE : 


The integer J specifies the material properties to be used 
for the inner n-leg thermoelectric washer segment (See 
Figure 1). The code numbers corresponding to each type 
of thermoelectric material whose properties are built into 
the code are as follows: 

J = 1; TEGS-3N material. 

J = 2; TEGS-2N material. 

J = 3; TEGS-3P material. 

J = 4; TEGS-2P material. 

J = 5; Ternary n-type material. 

J ~ 6; Ternary p-type material. 

The integer J specifies the materials properties to be used for 
the outer n-leg T/E washer segment (J defined as for NZ(1) 
above). 

The integer J specifies the materials properties to be used for 
the inner p-leg T/E washer segment (J defined as above). 

The intejer J specifies the materials properties to be used for 
the outer p-leg T/E washer segment of the P-leg (J defined as 
above). 

Output control parameter. Standard output format used 
except if: 

I RITE = I; NZ and Z array output suppressed. 

I RITE = 2; Radial temperature profile and temperature drop 
which are standard output for non-parametric calculations, 
are suppressed. 
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TABLE 1 (Continued) 


NZ(6) = NCLDH: 


NZ(7) = NCLDC: 


NZ(8) = KFTEMP: 


NZ(9) = IZ9: 


NZ(10) = NC: 


NZ(ll) = NGT1 : 


I RITE = 3; NZ and Z array; Radial temperature output 
suppressed. 

IRITE = 4; Output restricted to one page of input temperatures 
and calculated parameters per case. 

Parameter specifying inner clad material. 

NCLDH = 1; Stainless Steel 316 properties used. 

NCLDH = 2 or 0; Inconel 718 properties used. 

NCLDH = 3; Ta-lOW properties used. 

Parameter specifying outer clad material. 

NCLDC = 1 or 0; Stainless Steel 316 properties used. 

NCLDC = 2; Inconel 718 properties used. 

Parameter used to specify input and output temperature units. 
KFTEMP -= 0; Temperatures specified and listed in °K. 

KFTEMP f 0; Temperatures specified and listed in °F. 
Dimension input control parameter. 

IZ9 = 0; Module outer radius to be specified in Z(9) - See 
Table 2. 

IZ9 / 0; Radial thickness of outer T/E washer segment 
(r^ - r^, in Figure 1) to be specified in Z(9) - See Table 2. 
Number of axial sections (complete unit couples as shown 
in Figure 1) in the module (300 maximum). NC = 1 for 
all parametric studies in which module performance is based 
on results of a unit couple operating at average clad temper- 
ature conditions. 

Parameter specifying type of calculations to be done. 

NGT1 = 0; Open circuit and matched load calculations for 
a module in which all axial sections (see NZ(10), above) 
are connected in electrical series and thermal parallel. 
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TABLE 1 (Continued) 

NGT1 = 1; Open circuit calculations only for a module in 
which all axial sections (see N2(10), above) are connected 
in electrical series and thermal parallel. 

NGT1 = 2; Matched or fixed load calculations only for a 
module in which all axial sections (see NZ(10), above) are 
connected in electrical series and thermal parallel. 

NGT1 = 3; Open circuit calculations made for each 
individual couple (axial section) with no electrical or 
thermal connections between couples. 

NGT1 = 4; Matched load or fixed load (see 2(19) below) 
calculations made for each individual couple (axial section) 
with no electrical connections between couples. 

NGT1 = 5; Matched or fixed load (see 2(19) below) and 
open circuit calculations made for each individual couple 
as in NGT1 = 3 and NGT1 = 4 above. 

NGT1 = 6; Optimum load calculations (for maximum 
efficiency) made for each individual couple (axial section) 
with no electrical or thermal connections between couples. 

NGT1 = 7, Optimum load and matched load calculations 
made for each individual couple as in NGT1 = 6 and 
NGT1 = 4 above. 

N2(l 2) = PCMULT: Number of couples in the module for which individual couple 

calculations have been specified (N2(ll) > 3, above). 
Module performance is determined by multiplying appropriate 
parameters calculated for individual couples by PCMULT. 
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TABLE 1 (Continued) 


NZ(13) = INT: 


NZ(14) = INTTC: 


NZ(15) = iCRH: 


Results are strictly valid only for modules operating with 
uniform hot and cold clad temperatures. 

Method by which hot and cold clad temperatures are to be 
read as input. 

INT =• 0; Input hot and cold clad temperature pairs for each 
of the NC axial sections (see NZ(10) above). 

INT = 1; Axial hot and cold clad temperature profiles 
specified in previous case are used (for use in parametric 
studies). 

INT y 1; Input INT hot and cold clad thermocouple readings 
and interpolation will be performed based on axial locations 
(see NZ(14)) to determine hot and cold clad temperatures 
at the midpoints of each of the axial sections. This option 
can be used only if NZ(11)<(3. 

Parameter specifying method of reading axial locations of 
thermocouples (used only if NZ(13)>1). 

INTTC = 0 Input INT see (NZ(13)) hot and cold clad 
thermocouple readings. Individual thermocouples assumed 
to be uniformly spaced along circuit length of module with 
the first and last couples located at either end of the circuit. 
INTTC f 0 Input an axial distance corresponding to each 
thermocouple pair. Axial distances to be entered as specified 
in Section E below. 

Parameter specifying inner conductor ring material. 

ICRH = 1 Iron properties used. 

ICRH = 2; Tungsten properties used. 

ICRH f 1 and f 2: Molybdenum properties used. 
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NZ(33) = 1PUN: 


N2>(34) = IPDWT : 


N2(35) = ITHQ: 


NZ(36) = ITCR: 


TABLE 1 (Continued) 

Parameter specifying punched output in perturbation 
calculations. 

I PUN = 0; No punched output. 

IPUN = 1 ; A card of output parameters punched for each 
calculational case. Parameters listed are: Inner and outer 
diameters, average inner and outer clad temperatures, 
circuit length, load voltage, internal resistance, overall 
efficiency and power output in a IX, 2F7.3, 2F7.1, 5E10.3 
format. 

Parameter specifying power density or weight calculation 
output in perturbation routine calculations. 

IPDWT = 0; Power density (watt/cc) calculations printed. 
IPDWT 0; Weight calculations printed. 

Parameter specifying operating conditions for life test 
calculations (See Section IV). 

ITHQ = 0 or I; Fixed inner (hot) clad temperature. 

ITHQ = 2; Fixed or decaying heat input. 

Parameter specifying heat rejection conditions for life test 
calculations (See Section IV). 

ITCR = Oorl; Fixed outer (cold) clad temperature. 

ITCR = 2; Fixed radiator. 
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TABLE 1 (Continued) 


NZ(16) = ICRC Parameter specifying outer conductor ring material. 

ICRC = 1 Iron properties used. 

ICRC = 2; Tungsten properties used. 

ICRC f 1 and/ 2; Molybdenum properties used. 

N2(17) Not used. 

NZ(18) = IPIN: Parameter specifying material used for power lead pins (pins 

i 

extending through retainer rings at each end of module/ to 
which load circuit is connected). 

IPIN = 1; Iron properties used. 

IPIN = 2; Tungsten properties used. 

IPIN = 3; Molybdenum properties used. 

IPIN = 4or0; Nickel properties used. 

NZ(19), NZ(20), NZ(21) Not used. 

NZ(22)* = NODUMP: Control parameter used to request intermediate calculated 

parameters as output. 

NODUMP = 0; No intermediate information printed out. 
NODUMP = 1; A page of intermediate parameters (inter- 
face resistances. Joule heat, etc.) printed out after each pass 
through subroutine couple. 

NODUMP = 2; Intermediate parameters are printed out 
after the final pass through subroutine couple for each axial 
section. 

NZ(23) = MAXIND: If module current closure is not obtained after MAXIND 

iterations, calculation will be terminated. If zero is entered, 
MAXIND is set equal to 10. 
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TABLE 1 (Continued) 


N2(24) 


NZ(25)* 

NZ(26)* 

NZ(27) 


MAXTEM: if the criterion for temperature closure is not met for any 

couple after MAXTEM iterations, calculations will be 
terminated for this case and a dump of selected parameters 
will be given. If zero is entered, MAXTEM is set equal to 10. 

NPUMP: Parameter used to specify pump module parametric calculations 

(See section IV). 


NPERT: 

NR5J0PT 


Increment on number of couples to be used in parametric 
calculations (See Section IV). 

Parameter used to allow direct calculation of optimum T/E 
washer segmenting radius (r,. in Figure 1) in parametric 
calculations. 


NZ(28), NZ(29), NZ(30) Not used. 

NZ(31)* = NOPTIM Parameter used to specify temperature, load resistance, or 

geometry parametric calculations. (See Section IV). 

NOPTIM = 0 Perturbation subroutine (OPTIM) not entered. 
NOPTIM = 1 Standard value for perturbation calculations. 
One line of output listed for each set of temperatures, load 
resistance and dimensions. 

NOPTIM = 2 Used for T/E washer segmenting radius 
perturbations. Output listed only for optimum valve of r^. 
NOPTIM = 3 Temperature derivative option (See Section IV). 


NZ(32) = DRINCR Parameter specifying T/E washer segmenting radius pertur- 

bations. DRINCR corresponds to the increment to be applied 
to the inner segment thickness (in percent of the total 
specified T/E washer thickness). 
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FLOATING POINT INPUT PARAMETERS 


All floating point input, except for axial temperatures, are read into a list 
labeled Z. Table 2 lists the parameter to be read into each location and the method 
of entering data into this array is given in the section entitled DATA CARDS. The radi 
(r's) and axial dimensions referred to in the table are shown in Figure 1. 

All dimensions are in inches, resistances are in ohms and temperatures are in 
°K or °F depending on the value entered in NZ(8). Any exceptions are specified in 
Table 2. 

All entries designated by an asterisk in Table 2 refer to parameters which are 
zeroed at the beginning of each calculational case to avoid potentially expensive (in 
terms of computer time) errors. If non-zero entries are desired in change cases for any 
of these parameters* the values must be reset in each change case. All other entries 
in the Z array will maintain their previous values for all change case, unless changed 
by entering new values in subsequent cases. 
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TABLE 2 

TEMOD INPUT FLOATING POINT CONTROL PARAMETER DEFINITIONS 


Z0) 

= R(l) 

2(2) 

= DR(1) 

2(3) 

= DR (2) 

2(4) 

= DR(3) 

2(5) 

= DR(4) 

2(6) 

= DR(6) 

2(7) 

= DR (7) 

2(8) 

= DR(8) 

2(9) j 

= R(9) 


= DR(5) 

2(10) 

= ZLI 

2(11) 

= ZLN 

2(12) 

= ZLP 


2(13) = HH 

2(14) = HC 

2(15) = TOLTEM: 


Inside radius (r^) - inches. 

Radial thickness of inner clad (r^ - r^). See also 2(24). 

Radial thickness of inner insulating sleeve (r^ - r^). 

Radial thickness of inner conductor ring (r^ - r^). 

Radial thickness of inner T/E segment (r^ - r^). 

Radial thickness of outer conductor ring {r n - r ). 

/ o 

Radial thickness of outer insulating sleeve (r Q - r_). 

o / 

Radial thickness of outer clad (r Q - r_). See also 2(25), 

v o 

Outside radius r^, if NZ(9) = 0. 

Radial thickness of outer T/E segment, (r, - r c ), if NZ(9) / 0. 

6 O 

Axial length of insulating rings. 

Axial length of n-leg, T/E washer. 

Axial length of p-leg T/E washer. If zero is entered in this 
location, the axial length of the p-leg is calculated to 
optimize efficiency. See also 2(29). 

Thermal contact coefficients of hot insulating sleeves 
interfaces (Typically zero). 

Thermal contact coefficients of cold insulating sleeves 
interfaces (Typically zero). 

Tolerance to which temperatures must agree from one 
temperature iteration to the next to meet convergence 
requirement (TOLTEM = .001 if zero is entered) - See 
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2(16) = TOLCUR: 

2(17) = TPIN: 

2(18) = RPIN: 

2(19)* = RLOAD: 

2(20) = RECONC: 

2(21) = RECONH: 

2(22) = ZKEND: 

2(23) = RRREF: 


TABLE 2 (Continued) 

Equation 41 of WANL-TME-1 906. 

Tolerance to which current must agree from one current 
iteration to the next to meet convergence requirement. 
(TOLCUR = .001 if zero is entered) - See Equation 45 
of WANL-TME-1 906. 

Temperature at which electrical resistance of power lead 
pins (2(18)) has been evaluated. TPIN = 700° K (800° F) 
if 2(18) = 0. 

Electrical resistance of power lead pins connecting T/E 
circuit to external power leads. 

Load resistance. RLOAD = Internal generator resistance, 
i.e., matched load, if zero is entered and load calculations 
are specified. 

Contact resistivity coefficient of cold iron connector interfaces 
(See K mT in Equation 1 of WANL-TME-1 906). Typically 
RECONC = 0. 

Contact resistivity coefficient of hot iron connector interfaces 
(See ^1^-]- > n Equation 1 of WANL-TME-1 906). Typically 
RECONH = 0. 

Overall conductance of both module end closures (watts/°K 
if NZ(8) = 0; watts/ °F if NZ(8) / 0). 

Inner- to -outer clad radius ratio (r^/r^) at which end 
closure conductance, ZKEND, is evaluated. RRREF is used 
to scale 2KEND in parametric analyses in which radial 
dimensions are varied. 
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TABLE 2 (Continued) 


2(24) 

= CID: 

Ratio of inner clad thickness to module inner diameter. 

This parameter is used for parametric analyses in which the 
I.D. is being varied and is used to calculate DR(1) for each 
case. 2(2) must be set equal to zero. 

2(25) 

= COD: 

Ratio of outer clad thickness at module outer diameter. This 
parameter is used for parametric analyses in which the modul 
O.D. is being varied and is used to calculate DR(8) for each 
case. 2(8) must be set equal to zero. 

2(26) 

= ZNPERT: 

Increment for the n-leg axial thickness, 2LN, in parametric 
calculations (See Section IV). 

2(27) 

= itN/viiN: 

minimum n-ieg axial thickness, 2LN, in parametric calcula- 
tions (See Section IV). 

2(28)* 

= QREQ: 

Required module heat input rate (watts) in parametric 
calculations in which total heat input is specified. 

2(29) 

= 2LNLP 

Specified 2LN/2LP ratio, used only if 2(12), 2LP, is 
set equal to zero. 

2(30) 

= PEREQ: 

Specified power output for parametric calculations (See 
Section IV). 

2(31) 

= VREQ: 

Specified module load voltage for parametric calculations. 

2(32) 

= ZLREQ: 

Specified total module length for parametric calculations. 
(2(32) = 0 if couple axial dimensions are specified). 

2(33) 

= THINC: 

Increment for T u in temperature perturbation calculations. 

M 

2(34) 

= TC1NC: 

Increment for T- in temperature perturbation calculations. 
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TABLE 2 (Continued) 


2(35) = THMAX: 

2(36) = TCMAX: 

2(37) = DTMIN: 

2(38) = DRMIN: 

2(39) = RIINC: 

2(40) = R0INC: 

2(41) = R1 MAX : 

2(42) = R9MAX: 

2(43) = PWRDEN: 

2(44)* = QINPC: 


Maximum to be used in temperature perturbations. 
Maximum T^. to be used in temperature perturbations. 
Minimum radial temperature drop (T^ - T^,) for which 
calculations are to be performed. 

Minimum T/E washer radial thickness for which calculations 
are to be performed. 

Increment for module inner radius. 

If N2(9) = 0, ROINC is the perturbation increment 
applied to the module outer radius, r^. 

If N2(9) / 0, ROINC is the perturbation increment applied 
to the outer T/E washer radial thickness, AR^. 

Maximum module inner radius, r^, for which calculations 
are to be performed. 

If N2(9) = 0, R9MAX is the maximum module outer radius, 
r^, for which calculations are to be performed. 

If N2(9) / 0, R9MAX is the maximum outer T/E washer 
radial thickness, AR^, for which calculations are to be per- 
formed. 

Fuel specific power density (watts/cc). Fuel is assumed to 

fill volume enclosed by the I.D. of the module inner clad. 

T u will be calculated for cases in which PWRDEN > 0. 
tl 

Specified heat input rate (watts) for each couple. T|_| will 
be calculated for cases in which QINPC > 0. 
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2(45) = QGAM: 


2(46) 

RLPERT: 

2(47)* 

RLMAX: 

2(48), 2(49) 


2(50)* 

2LIBRH 

2(51 )* 

2LIBRC 

2(52), 2(55) 


2(56) 

HRINC: 

2(57) 

HRMAX: 

2(58) 

HFLF: 

2(59) 

TREJ: 


2(60) 


TABLE 2 (Continued) 

Amount of gamma heat generated within the lead telluride 
of each couple in an internally fueled gamma emitting isotope 
configuration. This parameter must be determined by an 
independent gamma heat shielding analysis. 

Increment for load resistance (2(19) applied in perturbation 
subroutine. No perturbation performed if RLPERT = 0. 

Maximum load resistance for which calculations are to be 
performed in load resistance perturbation subroutine. 

Not used. 

Axial thickness of the tungsten foil diffusion barriers incorpor- 
ated in the insulating washers in annulus D of Figure 1. 

Axial thickness of the tungsten foil diffusion barriers incorpor- 
ated in the insulating washers in annulus E of Figure 1. 

Not used. 

Increment (hours) to be used in time increment calculations 
(See Section IV). 

Maximum time (hours) to be used in time increment calculations. 

Half life (years) of isotope fuel in time increment calculations. 

Ambient heat rejection temperature to be used in time 
increment calculations performed for fixed radiator configur- 
ation (See Section IV). 

Not used. 
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2(61) = WTCON : 

Ratio of module end closure weight to cross sectional area 

.2 

of module (pounds per in ). 

2(62), 2(70) 

Not used. 

2(70 + J)** 

Percent increase desired in Seebeck coefficient of Jth 
thermoelectric material where J = 1 to 6 (See N2(l)). 

2(80 + J)** 

Percent increase desired in resistivity of the Jth thermoelectric 
material where J = 1 to 6 (See N2(l)). 

2(90 = J)** 

Percent increase desired in thermal conductivity of the Jth 
thermoelectric material where J = 1 to 6 (See NZ(1)). 


** 


These entries refer to percent increase per 1000 hours in life test calculations to 
account for module degradation as a function of time as discussed in Section IV 
of this report. 
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E. DATA CARDS, ORDER AND FORMATS 

1 . Comment Card 

The first data card read into the program is a comment card. This comment 
should identify the case being run and will print out at the top of each page of output. 

The card is read using a 18A4 format so that the comment is restricted to a maximum of 
72 spaces and should be centered in the 72 space field on the data card. 

2. Fixed Point Data for NZ Array 

Immediately after the comment card? fixed point data are read. A maximum 
of eleven (1 1) entries in the NZ (fixed point data)list may be made on one card. The 
format used in reading fixed point data is FORMAT (312, 11 (3X, 13)). 

The first two digit integer refers to the number of pieces of data given 
on the card. This must be a number form 0 to 1 1. Any entry other than zero in card columns 
3-4 indicates that no more fixed point data will follow the card being processed. The 
two digit integer in card columns 5-6 indicates the NZ subscript of the first piece of 
data on the card being processed. The remaining input data centered on each card are 
entered sequentially into the NZ array. 

3. Floating Point Data for Z Array 

After processing a fixed point data card with an entry in card columns 3-4, 
the floating point data are read. A maximum of six (6) entries in the Z list (floating 
point data) may be made on one card. The format used to read this data is: FORMAT 
(312, 6F12.6). The first two digit integer refers to the number of pieces of data given 
on the card. This must be a number from 0 to 6. Any entry other than zero in card columns 
3-4 indicates that no more floating point data is to be entered into the Z list after the 
present card is processed. The two digit integer in card columns 5-6 indicates the Z 
subscript of the first piece of data on the card being processed. The remaining data on 
the card are entered sequentially into the Z array. 
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4. Temperature Data 

The method used to input clad temperature data is determined by the control 
variable entered in NZ(13). All temperatures must be entered in Kelvin units. Three 
options are available for specifying temperature data. 

a. NZ(I3) = 0. If NZ(13) has been set equal toQ, hot and cold clad tempera- 
ture pairs for each axial section must appear immediately after the first 
floating point data card with a non-zero entry in card columns 3-4. 

The number of temperature pairs to be read is the number previously 
entered in NZ(10). 

Six temperatures pairs are read per card using a 12F6.5 Format. The 
first temperature of each pair must be the hot clad temperature. In the 
event that a zero is encountered in any of this temperature input data, 
an error message is printed and all calculations for the case are suppressed. 
Hence it is possible to make an input error for one case without affecting 
the following cases. 

b. NZ(13) =1. In the event that the effects of geometry variations are 
being investigated (as in parametric studies), it is often desired to run 
many cases using the same clad temperatures. If the number entered 
in NZ(13) is 1, no further data cards are read after the final card 
entering data in the Z array. The code uses the temperature data entered 
in the previous case. 

c. NZ(13) y 1. This option is used to calculate the performance of a module 
for which experimental clad temperature measurements have been made 

at various locations along each clad. If NZ(13) has been assigned greater 
than 1, the computer will interpolate on input hot and cold clad thermo- 
couple readings to determine the hot and cold clad temperatures of each 
axial section. The number entered in NZ(13) corresponds to the number 
of hot and cold clad thermocouple pairs (not to be confused with the 
number of thermoelectric couples in the module) to be read. 
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In order to perform the interpolation, of course, it is necessary to specify 
the axial location of each of the thermocouples. These locations are 
specified as the distance (in inches) from the leading edge of the first 
axial section of the module. 

Since standard thermocouple instrumentation is used for most modules 
tested in the Compact Thermoelectric Converter program, the axial 
locations of these standard thermocouples are built into the code. By 
assigning NZ(14) = 0, the thermocouples are assumed to be uniformly 
spaced along the circuit length of the module with the first and last 
couples located at either end of the circuit. In this case, the first card 
(or set of cards) after the Z array data cards must contain the hot clad 
temperature data. The next card (or set of cards) must contain the cold 
clad temperature data. All of these cards use a 12F6.5 Format. 

To specify externally the thermocouple locations, the number of thermo- 
couple pairs must be the number entered in NZ (13). If NZ(14) has 
been assigned any value other than zero (and if NZ(13) has been assigned 
a number greater than 1) the computer will begin to read axial locations 
for each thermocouple pair immediately after reading the last floating 
point data card. 

After reading the thermocouple locations, the computer will read I NT 
(the number entered in NZ(13)) hot clad thermocouple temperatures and 
the INT cold clad temperatures. Each of the three lists begins on a new 
data card and each card uses a 12F6.5 Format. 

A zero entry in either the hot or cold clad thermocouple data is used to 
indicate the absence of a thermocouple reading at a particular location. 
When performing the interpolation to determine the temperature at the 
mid-point of each axial section, all zero entries in the temperature lists 
are ignored. 
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5. Multiple Case Runs 

Multiple runs can be made using TEMOD by simply stacking sets of input 
data. After reading all the input data for a given case, the calculations are performed 
and the output listed. The computer returns to input to search for another comment card. 

If none exists, calculations are terminated. 

With the exceptions noted in Sections C and D above, entries in the Z array, 
NZ array and axial temperature profiles are not zeroed from one case to the next. Hence, 
it is normally not necessary to reread any data which has not changed from one case to 
the next. Only the parameters which are changed from the previous case need be read 
in. 

Care must be taken, however, to include a comment card, at least one 
card with fixed point data and one with floating point data for each case. If, for example, 
no floating point parameters change from one case to the next, a dummy card with only 
an entry in card column 4 must be included in the data cards in place of the floating point 
data card. 
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IV. PARAMETRIC CALCULATIONS 

A. GENERAL 

The TEMOD code has been written to allow various type of parametric calculations 
to be performed. Incrementation of parameters (i.e. operating temperatures, load resistances 
and component geometries) is done automatically by the code eliminating the need for 
voluminous input data decks. In addition, the output for these types of calculations is 
restricted to one line per case. Provisions have been made to allow specifications of 
module operating parameters such as electrical power, voltage, and/or heat input to 
meet required module design operating conditions. These conditions are met internally 
using program logic to calculate the required number of thermoelectric couples and/or 
module circuit length. 

In each case, the calculations are made on a unit couple (see Figure 1). Thus, 
in each parametric calculation case, the control parameter NZ(10) should be set equal 
to 1 and one set of temperature pairs should be entered. 

The control parameters, geometry and operating temperatures for the first 
case of each parametric calculation must be set in accordance with instructions given in 
Section III of this report. The pertinent control parameters which must also be set for 
each type of parametric calculation is discussed below. 

B. TEMPERATURE PARAMETRIC WITH SPECIFIED GEOMETRY 

This option is normally used to determined performance of a specific type of 
module operating over a wide range of average hot and cold clad temperatures. In addition 
to the control parameters required to specify the module component materials and dimensions* 
the following entries are required to perform this type of parametric: 

1. Set NZ(31) = 1. 

2. Set Z(31) = Z(32) = 0. 

3. Set Z(33) through Z(36) equal to the appropriate values as listed in Table 2. 

4. Set the initial temperature pair at the lowest hot and cold clad levels of 
interest in the parametric. 
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c. 


LOAD RESISTANCE PARAMETRIC WITH SPECIFIED MODULE 


This load is used to calculate load curve characteristics of a specific type of 
module at a specific set of operating temperatures. The result?, of course, will show that 
the power output will approach zero as the load resistance approaches either limit (zero 
or infinity) and that maximum power occurs at the point where the load resistance very 
nearly equals the generator resistance, i.e. matched load. (For calculations performed 
at fixed hot and cold junction temperatures, maximum power would occur precisely at 
matched load. However since the module performance calculations are made at fixed 
clad temperatures and since the junction temperatures do vary as the load resistance varies, 
the maximum power point does not occur precisely at matched load. Since the deviation 
between matched load power and maximum power is extremely small, no distinction is 
made between the two and an option has been built into the code to allow a direct 
determination of the matched load performance eliminating the need for running a load 
parametric.) 

The load parametric calculations will also indicate that the load resistance 
at which module efficiency is optimized (defined as "optimum load") is approximately 
20 to 30 percent higher than the module internal resistance. An option has also been 
built into the code to allow a direct determination of optimum load performance without 
running a load parametric. 

Load resistance parametrics may be specified over a narrow or a wide range 
of resistance values. A parametric with a step size and increment range resulting in 
more than 200 separate cases is considered a wide range parametric. In the wide range 
calculations, the increment size is increased by a factor of ten after every ten calcula- 
tions. This allows the parametric to cover an extremely wide range of resistance with 
only ten calculations for each order of magnitude. 

In addition to the control parameters and input data required to specify the 
module component materials and dimensions, the following entries are required to perform 
the two types of load resistance parametrics: 
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1. Narrow range parametric: 

a. Set NZ(31) = 1 

b. Set Z(19) equal to the initial (lowest) load resistance level. 

c. Set Z(46) and Z(47) equal to the appropriate values. 

2. Wide range parametric: set all parameters as described in (1) above except 
Z(46) = 0. 

D. TEMPERATURE PARAMETRIC WITH SPECIFIED VOLTAGE AND CIRCUIT LENGTH 
This option allows parametric calculations to be performed to determine optimum 

module dimension for applications in which a load voltage and total circuit length has 
been prescribed. For these calculations/ the number of thermoelectric couples required 
to meet the voltage specification at either matched or optimum load conditions is determined 
internally for each set of operating temperatures. 

In addition to the control parameters and input data required to specify the module 
component material and dimensions, the following entries are required to perform this type 
of parametric: 

1 . Set NZ(1 1) = 4 or 6 (matched or optimum load calculations). 

2. Set NZ(31) = 1 

3. Set Z(12) = 0. to allow p-leg axial thickness to be optimized with respect 
to n-leg thickness for each set of operating temperatures. 

4. Set Z(19) = 0. (A specification of load resistance, length, and voltage 
amounts to an overspecification of the module performance). 

5. Set Z(3?) through Z(36) equal to the appropriate values as listed in Table 2. 

E. PARAMETRIC ON NUMBER OF COUPLES WITH SPECIFIED LOAD VOLTAGE 
AND CIRCUIT LENGTH 

This type of parametric is very similar to D above, except the load voltage 
specification is met in each case by setting the load resistance equal to the appropriate 
value. This, of course, is done internally for each case by the program. This option is 
intended for use in applications where operating temperatures have been established by 
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system operating constraints, hence, temperature parametrics should not be attempted. 

1. Set NZ(12) = 1 (the minimum number of couples required to achieve the 
specified lead voltage is determined by program logic). 

2. Set NZ(26) = 1 (the normal desired increment on the number of couples). 

3. Set NZ(31) = 1. 

4. Set Z(12) = 0. To allow p-leg axial thickness to be optimized with respect 
to n-leg thickness for each case. 

5. Set Z (27) equal to the minimum n-leg thickness to be considered (this 
corresponds to a specification of the maximum number of couples and 
determines when the parametric calculations will be terminated). Z(2 7) = 
.020 inch if no value is input. 

6. Set Z(31) and Z(32) equal to the appropriate values as specified in Table 2. 

7. Set Z(33) through Z(36) equal to zero. 

F. PARAMETRIC ON N-LEG/P-LEG AXIAL THICKNESS RATIO WITH SPECIFIED 
CIRCUIT LENGTH 

The code has been written to allow a direct determination of the optimum n-leg/ 
p-leg axial thickness ratio. This option, however, allows a determination of the effects of 
varying this ratio on module performance. 

In addition to the control parameters and input data required to specify the 
module component materials and initial case dimensions, the following entries are required 
to perform this type of parametric: 

1. Set NZ(12) equal to the number of couples in the module. 

2. Set NZ(26) equal to the number of washer thickness perturbations required. 

3. Set NZ(31) = 1 . 

4. Set Z(11) and Z(12) equal to the initial n- and p-leg washer axial thicknesses 

5. Set Z (25) equal to the length by which the n-leg washer thickness is to 

be increased and the p-leg washer thickness decreased in each perturbation. 


D-28 




Astronuclear 

Laboratory 


G. 


PARAMETRIC ON NUMBER OF COUPLES WITH SPECIFIED LOAD VOLTAGE 


HEAT INPUT OR POWER OUTPUT 

In many module applications/ system constraints govern the module voltage and 
heat input or power output requirements. For a given set of operating temperatures, 
the optimum module meeting these requirements can be determined using a parametric 
routine built into the code. In this routine, the number of couples in the module is 
increased from the minimum number required to produce an open circuit voltage greater 
than the required voltage to a maximum number which is determined primarily by 
fabrication limits (minimum thermoelectric washer thickness). The total thermoelectric 
circuit length is dictated, primarily, by the heat input or power output specification. 

In addition to the control parameters and input data required to specify the 
module component materials and initial case dimensions, the following entries are required 
to perform this type of parametric: 

1. Set NZ(12) = 1 (the minimum number of couples required to achieve the 
specified load voltage is determined by program logic). 

2. Set NZ(26) = 1 (the normal desired increment on the number of couples). 

3. Set NZ(31) = 1. 

4. Set Z(12) = 0 to allow p-leg axial thickness to be optimized with respect 
to n-leg thickness for each case. 

Set Z(28) equal to the desired module thermal power input, or 
Set Z(30) equal to the required electrical power output. 

6. Set Z(31) equal to the desired load voltage. 

7. Set Z(32) through Z(36) equal to zero. 


1 : 


H. 


PARAMETRIC TO DETERMINE OPTIMUM SEGMENTING RADIUS, 


For many applications, it is desirable to use radially segmented thermoelectric 
washers in either the p- or n-legs of the module. Thermoelectric materials can be doped 
at different levels to provide maximum conversion efficiency in different operating 
temperature ranges. Using a material doped for optimum performance at high temperatures 
in the inner segment of each thermoelectric washer, and an alternate material composition 
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providing optimum performance at lower temperatures in the outer segment of each thermo- 
electric washers can provide efficiency improvements. Given two types of n-type and/or 
two types of p-type material, then, there is an optimum radius at which the materials can 
be segmented. An option has been provided to allow a determination of this optimum 
segmenting radius. 

Since the optimum segmenting radius is a function of the module operating tempera- 
tures, this parametric option may be used in conjunction with the temperature parametric 
option discussed in (B) above. 

In addition to the control parameters and input data required to specify the module 
component materials and initial case dimensions, the following entries are required to perform 
this type of parametric: 

1 . Set NZ(31 ) = 0. 

2. Set NZ(32) equal to the increment to be applied to the inner segment 
thickness in percent of the total specified T/E washer thickness. Normally 
NZ(32) = 5 provides adequate caiculational resolution to determine the optimum 
segment radius. 

3. Set 2(5) = 0. This initializes the segmenting radius such that the radial 
thickness of the inner T/E washer segment is zero. 

4. Set 2(28) through 2(32) equal to zero. 

5. Set 2(33) through 2(36) equal to the appropriate values listed in Table 2 
to allow operating temperature perturbations. 

6. Set the initial temperature pair at the lowest hot and cold clad levels of 
interest in the parametric. 

With this deck setup, a line of output will be listed for each segmenting radius. 

The first case (or line of output) for each temperature pair will correspond to a zero thickness 
outer T/E washer segment, i.e. 

l °°( r 6 " r 5 ) 

= 0 pet. (See Figure 1). 

r 6 ‘ r 4 
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or 


100 (r - r ) 

- 1 =100 pet. 

r 6 " r 4 


The final case for each temperature pair will correspond to a zero thickness 
inner T/E washer segment. 


100 ^ 

r 6 ' r 4 


= 100 pet. 


or 


1 00 (r 5 - r 4 ) 
r 6 " r 4 


0 pet. 


Since normally interest is restricted to the optimum segmenting radius only, 
an alternate form of this parametric can be used to reduce both the amount of output and 
the required computer time to perform the parametric. This option is specified in a manner 
identical to that discussed above except NZ(31) is set equal to 2 . With this option 
incrementing of the inner washer thickness is performed until a maximum efficiency has 
been obtained, a line of calculated parameters corresponding to this optimum geometry is 
printedf and the temperatures are then incremented and the procedures repeated. 


I. TEMPERATURE DERIVATIVE CALCULATIONS 

During the course of reduction and analysis of experimental data from modules 
being tested, it has been found necessary to eliminate performance variations resulting 
from minor temperature fluctuations. For each experimental data set, it is desirable to 
determine analytically what the performance parameters would have been if the module 
had been operated at the exact prescribed operating temperatures. 

This determination can be made by determining the derivatives of the primary 
performance parameters (i.e. effective Seebeck coefficient, a, internal resistance, Rg, 
and thermal impedance Tl) with respect to hot and cold clad temperatures. 
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A routine has been provided in the code to allow this evaluation. Calculations 
are performed at 25°F increments on either side of the design hot and cold T^, 
clad temperatures. The parameters 



are computed and listed on the basis of these calculations. In these expressions, the effective 
Seebeck coefficient is defined: 


V 

oc 



where V c is the module open circuit voltage, and the thermal impedance is defined 


Tl = 


t h' t c 

Q 


where Q is the total module heat input. 

In addition to the control parameters and input data required to specify the module 
component materials and initial case dimensions the following entries are required to 
perform this type of parametric: 

1. Set NZ(ll) = 4. (The calculations are typically performed under matched 
load conditions except in the case of open circuit tests.) 

2. Set N2(31 ) = 3. 

3. Set 2(28) through 2(36) = 0. 

4. Set the input temperatures at the design operating levels of the module. 

J. RADIAL GEOMETRY PARAMETRICS 

Quite often it is desired to perform any of the parametric options discussed above 
for a family of modules with varying overall radial geometries. For this reason an option 
is provided in the code to allow the radial geometry incrementation to be handled internally 
by the code. 
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Provisions have been made to allow perturbation of the module inner radius and 
module outer radius or lead telluride radial thickness. Both of these geometry incrementing 
options may be used in conjunction with any of the parametric routines discussed previously. 

1 . Inner Radius/Outer Radius Parameters 

In addition to the control parameters and input data required to specify the 
module component materials and initial case dimensions, along with the entries to perform 
any of the previously discussed parametric calculations, the following entries are required 
to perform the radial geometry parametrics: 

a. Set N2(9) = 0 (See Table 1). 

b. Set 2(1) equal to the smallest inner radius of interest. 

c. Set 2(2) = 0. Inner clad thickness should be scaled linearly with module 
inner radius as specified by 2(24). 

d. Set 2(8) = 0. Outer clad thickness should be scaled linearly with module 
outer radius as specified by 2(25). 

e. Set 2(9) equal to smallest outer radius of interest. 

f. Set 2(23) through 2(25) equal to the appropriate values as specified in 
Table 2. 

g. Set 2(38) equal to the minimum T/E washer radial thickness for which 
calculations are to be performed. This is necessary since there may be 
rj, r^ combinations in the range of interest corresponding to negative 
or very small T/E washer thicknesses. A substantial amount of computer 
time may be required to achieve temperature convergence on these cases 
of little or no interest. 

h. Set 2(39) through 2(42) equal to the appropriate values as specified in 
Table 2. 

i. Set 2(61) equal to the appropriate value as specified in Table 2. 
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2. Inner Radius/Thermoelectric Washer Radial Thickness Parametrics 

This option is very similar to (1) above except that the outer radius of the 
module is controlled by the specified T/E washer thickness. All parameters should be set 
as discussed in (1) above except the following: 

a. Set N2(9) = 1 (See Table 1). 

b. Set 2(9) equal to the smallest T/E washer radial thickness of interest. 

c. Set 2(38) = 0. 

d. Set 2(42) equal to the maximum T/E washer radial thickness of interest. 

K. PUMP MODULE PARAMETRIC STUDIES 

An option has been built into the TEMOD code to allow a determination optimum 
dimensions for tubular modules designed to provide electrical power for electromagnetic 
pumps. Since electromagnetic pumps require high current/low voltage power, typical 
pump modules have extremely low internal resistances. For this reason, inner and outer 
conductor ring radial thicknesses are very critical and optimum thicknesses must be 
determined parametrically. 

The parametric routine included in the code requires a specification of required 
current and load voltage along with operating temperatures. The radial thicknesses of 
the inner and outer conductors and thermoelectric washers are each varied independently. 
Axial dimensions for each case are calculated on the basis of meeting the required current 
and voltage specifications. 

In addition to the control parameters and input data required to specify the module 
component materials and initial case dimensions, the following entries are made to perform 
this type of parametric: 

1. Set N2(9) = 0. (See Table 1). 

2. Set N2(25) equal to the desired number of perturbations to be performed 
on each conductor and T/E washer radial thickness. 

3. Set 2(4), 2(5), 2(6) and 2(9) equal to the appropriate minimum values of 
interest. 
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4. Set 2(8) = 0. Since module outer radius will vary in parametric, outer clad 
thickness should be scaled linearly with r^ as specified by 2(25). 

5. Set 2(12) = 0. Program logic will determine the axial thickness of the T/E 
washers required to achieve voltage and current specifications for each case. 

6. Set 2(19) = 0. Although the input voltage and current specify a load resistance, 
this value is calculated internally. 

7. Set 2(23) and 2(25) equal to the appropriate values as specified in Table 2. 

8; Set 2(30) equal to the specified power output (the product of the specified 

current and voltage). 

9. Set 2(31) equal to the specified voltage. 

10. Set 2(40) equal to the radial increment to be applied to the conductor rings. 

The T/E washer increment is half as large as the conductor ring increment. 

11. Set 2(61) equal to the appropriate value as specified in Table 2. 

L. CALCULATION OF PERFORMANCE AS A FUNCTION OF TIME 

During operation of thermoelectric generators various factors can produce performance 
variations as a function of time. Obviously, if the operating temperatures vary, the perfor- 
mance will be affected. These temperature variations can be caused by externally controlled 
modification or by a decay of the fuel in the case of an RTG application. In the latter 
case, the module cold clad temperature will vary as the amount of heat to be radiated 
is reduced. 

In addition to operating condition variations, module performance is affected by 
a degradation process. The effects of this process can be simulated by modifying the thermo- 
electric material properties in the appropriate manner. In lead telluride generators these 
effects are simulated in increasing the resistivity of the n-type thermoelectric material 
to compensate for the diffusion of tellurium into the material. The rate at which the resisti- 
vity increases is a function of operating temperatures and washer axial thickness, and must 
be given as input to the calculations. Provisions are also included for modifying any of the 
other themoelectric properties in a similar manner. 
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In addition to the control parameters and input data required to specify the module 
component materials and dimensions, the following entries are required to perform this 
type of calculation; 

1. Set N2(35) equal to the appropriate value discussed in Table I to specify 
either constant T|_| or heat input conditions. 

2. Set N2(36) equal to the appropriate value discussed in Table 1 to specify 
either constant T^. or fixed radiator calculations. 

3. Set 2(19) if fixed load resistance calculations are required (2(19) = 0 for 
matched load). 

4. Set 2(31) if fixed load voltage calculations are required. 

5. Set 2(56) equal to the desired time increment (in hours). 

6. Set 2(57) equal to the maximum time (hours). 

7. Set 2(58) equal to the isotope half life (in years) if heat decay calculations 

are desired. 

8. Set 2(59) equal to the heat sink temperature if fixed radiator calculations 
are desired. 

9. Set 2(71) through 2(100) equal to the appropriate T/E material property 
rate of change as specified in Table 2. Rates of change to specified as 
percent change per 1000 hours. 

10. Set temperatures at the beginning-of-life levels. 
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V. TEMOD FORTRAN PROGRAM LISTING 


PROGRAM TEMOD ( INPUT , OUTPUT , PUNCH# TAPE 5slNPUT,TAPE b-OUTPUT, 

1 TAPE 7*PUNCH) 

C TEMOD CHARLES ROSE.. ...BLDG. 2 

DIMENSION NZ(50) ,Z(100) ,XTHERM(30) , THERMc (30) , THERMH (30) , HD (b, 8 ) , 

1 HD l (8) ,HD2(8) ,HD3(lb,2) , ASTARl (2) ,ZOPMAX (5) , ZOPT (5) , ZOPT 1 (5) , 

2 CKF 1 (2) ,CKF2(2), TMKF (2) 

COMMON /T I TLE/FETM (3),CASCAD(8,2) , SSIN (3) , CMT (25)# COND (3, 5) , TEMKF , 
1 CRM0(2,2),I3N,I2N,I3P,I2P 
COMMON /MOLIF/HRS,HRINC.HRMAX,hFLF,RADPCT, ITCR.ITHQ 
COMMQN/TERIT/DTAV,OHM,QCV,RLOAD#QMOD, TCAV 

COMMON /MDCP0P/RAD(9) ,DR(8) ,TEmP(9,300) , V0C(300) ,RPC(300) ,PE(300) , 

1 QT (300) ,NOPTIM, VPC(300) ,RLPC(300) ,DTH0D(30O) ,CUR(300) ,LNRAD(9) , 

2 ZLP,ZLN,ZLI#ZLNLP,ZLPNI,ZLTE#nGT1,NGT2, ITPERT.PCMULT, NODUMP, I VO, 

3 QGEN,QGENL,GTE,U0UT,ZID,Z0D,ZkEND,ZK9,ZkR,RPPC,RPN,T0lTEM.TREJ» 

9 vreq, pereq.zlreq, wtcon.zlpn, IdUMZ , RN, DGRF , DGTRF , DGLRF , RADK , JDUMP , 
5 DP (30) , NC 

COMMON /M0DCPL/CURM0D(15), IT(300),DTTE,C0NRN,C0NRP,ZLlP,ZNC» 

1 ALRb4,MAXTEM,QGAM,CGAM, INDI,RtNT4,RINT6,RC0N4,RC0Nb,ClFEH,ClFEC, 

2 CONST 1.C0NST2.C0NST3. CONS r4,C0N3T5,NR50pT,SEGTE,ICRMDL,R ADC 
COMMQN/TEFO/NCLDC.NCLDH, I CRH, icRc, ipin,zlibrh,zlibrc,tpin.rpin 
COMMON /MODOPT / IPDWT , DR4 , DR I NcR , I R I TE , DRTE , JOPT I M , I Z9 , IENDK , I PUN , 

1 RIINC,ROINC,RIMAX,R9MAX,OTMIN,OREQ, THINC, TCINC, TCMAX.THMAX 
2,NPERT,ZNPERT,CKFP,CKFT,CTE0C,CTElC,DELTLf RLPERT,RLMAX,ZNMIN,Z9 SaV 
C0MM0N/TEPMP/NPUMP.DR8,C002 
EQUIVALENCE (H03(20) , ASTaRI (U) 

DATA HD/4HINNE, 4HR CL,4HAD , 4H , 4H , 4H , 

1 4HINNE,4HR IN,4HSULA,4hTING,4H RIN.4HG , 

2 4HH0T » 4h SID,4HE CO , 4hNDUC , 4HT UR , 4HR ING , 

3 4H ,4HINNE,4HR T/,4hE W A , 4HSHER , 4HS , 

4 4H , 4H0UTE , 4HR T/,4hE W A , 4HSHER , 4HS , 

5 4HC0LD , 4h SID, 4HE CU, 4hNDUC, 4HT0R , 4HRING, 

b 4H0UTE , 4HR IN , 4HSUL A , 4 hT I NG , 4H RIN.4HG , 

7 4H0UTE , 4HR CL,4HAD ,Ih ,1H ,iH / , TMKF / 1 HK , 1 HF / , K ASE 1 / 1 / , 

8 HD1 (2)/2HBN/,HDl (7)/2HBN/,HD2/3*lH ,2*lH*,3*lH /,HD3/20*lH ,IH*, 

9 b* 1 H ,2*1H*,3*1H /,NZ/50*0/,Z/100*0./,CKFl/0.,459.b/,cKF2/l,,l,8/ 
REAL LNRAD 

FK2K(T)«(T+CKFP)/CKFT 
CID»CL0CK(CMT(21),CMT(22) ,CMT(23)) 

93 READ (5,94) (CMT(K) ,K»3,20) , ASTOP 

94 FORMAT (18A4.7X.AI) 

Z(l9)*o. 

Z ( 2 8 ) * 0 . 

Z ( 44 ) «o » 
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Z (45) *0« 

Z(47) =0, 

Z(50)=0. 

Z(51)=0. 

Z(62)*0. 

. Z(67)=0. 

DO 90 1=70,100 

90 Z(IJsO. 

UGEN=0, 

ZKK=0, 

CURMOD(1)*0. 

NZ(20)=0 
nZ(22)sO 
N Z ( 25 ) s 0 
NZ(2fe)sO 
NZ (31) =0 
lTPERTsO 
I ND I = 0 
LNDUMPsO 
1VD=1 

IF CASTOR, NE. H02(l )) STOP 

c - read fixed point input data 

95 READ ( 5 » 9 1 ) N,L,J, ( IT ( I ) , 1*1 , N) 

91 FORMAT (312, 11(5 X, 13)) 

JsJ-1 

00 92 1=1, N 

92 NZ(ItJ)»IT (I) 

C - READ FLOATING POINT INPUT DATA 

IF (L.FQ.O.AND.N.NE.O) GO TO 95 

3 READ (5,2) N,L,J, ( VPC ( I ) , 1=1, N) 

2 FORMAT (IX, II ,2I2,6E12.6) 

J=J-1 

DO 4 I «1,N 

4 Z ( 1 + J) sVPC ( I ) 

IF (L.EQ.O.AND.N.NE.O) GO TO 3 

c - determine type of material specified for each t/e segment 

I3N*NZ ( 1 ) 

I2N = NZ ( 2 ) 

I3P*NZ(3) 

I2P = nZ ( 4 ) 

IRITE«nZ(5) 

NCLDHsnZ(S) 

NCLDCsnZ (7) 

KFTEM=1 
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IF (NZ(8),NE.O) KFTEM»2 
CKFPsCKFI (KFTEM) 

CKFT*CkF2(KFTEM) 

TEMKFsTMKF(KFTEM) 

1Z9« 1 

IF CNZ(9),NE,0) I Z9*2 
NCsMAX0(NZ(10),1) 

ZNC ■ NC 
NGTl«NZ(ll) 

PCMULTaAMAXO(NZ(12),l) 

IF (KASE1.EQ.1.AND.NZU3),EQ,1) NZ(13)»0 
INT * NZ ( 1 3) 

INTTCeN2(14) 

IC«H«NZ(15) 

ICHC«NZC16) 

IPIN«MlNQ(NZ(18),4) 

IF (IPIN.LE.O) I PI N«4 

NXTRP « NIN0(MAX0(NZ(19) , 1) ,5J 

SEGTE=10. 

IF (NZ(21).NE,0) SEGTE*NZ (21 ) 

NODUMPaNZ (22) 

JOUMPsnODUMP-1 
MAXINDsHAX0(NZ(23) ,10) 

MAXTEMsO 

IF CNZC2«),GT,10) MAXTEMaNZ(24)-lO 
NPUMP«nZ(25) 

NPER T a NZ (2b) 

NR5QPT aNZ (27 ) 

NOP T IMsNZ (31 ) 

IF (NOPTIM.GT.O) NC«1 
DRINCRa,01*FL0AT(NZC32)) 

XPUNaNZ (33) 

IF (NZ(34) ,NE.O) IPDWT *2 
ITHQaMAX0(NZ(35),t) 

ITCR«MAX0(NZ(3bJ, 1) 

ICRMOL*NZ (50) 

32 JOPTlHaO 


ISECT 

■1 


IPDWT 

«1 


RAO ( 1 

)• 

Z(t) 

OR 

(1) 

«Z 

(2) 

OR 

(2) 

*Z 

(3) 

OR 

(3) 

■Z 

(4) 

DR 

(4) 

«Z 

(5) 
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DR(6)*Z(6) 

DR(7)*Z(7) 

uwa = z c s j 

Z9SAVsZ(R) 

ZLlsAMAXl ( Z ( 1 0 ) , 1.0E-10) 

ZLN»Z(ll) 

ZLPsZCia) 

HH s Z ( 1 3 ) 

HC*2(U) 

TOLTEMsAMAXl (Z(15), ,001) 

IF (NGT1.EG.2) TOLTEMslO,*TOLTEM 
rOLCURaZ(lfe) 

IF (TOLCUR.LE.O.) TOLCURs.OOl 
TPIN8FK2K (Z (17) ) 

IF (TPlN.LE.256,) TPINb700, 

RPIN*Z(16) 

KLOADsZ (19) 

RECONCsZ (20) 

RECONHsZ (21 ) 

ENDZKsZ(22)*CKFT 
WRREF*Z (23) 

CI02*2,0*Z(2«) 

C002«2,0*Z (2b) 

ZNPERTsAMAXl (Z (26),. 001) 

ZNMINSAMAX1 (2(27) , .030) 

QREQsZ (28) 

ZLTE sZLP+ZLI+ZLl+ZLN 
IDUMZsl 

IF (ZIP.NE.O.) GO TO IS 
ZLTE sO, 

IDUMZsa 

ZLNLPsZ(29) 

IF (ZLNLP.GT.O.) IDUHZs3 
15 PEREQ*Z(30) 

VREQsZ (31 ) 

ZLREQsZ (32) 

THINCsZ( 33)/CKF1 
TClNCsZ(3i»)/CKFT 
THHAXsl 100, 

IF (THINC.GT.O.) T HM AX sF K 2K ( Z ( 35 ) ) +0 . 5* T H I NC 
TCHAXsO. 

IF (TCINC.GT.O.) TCMAX«FK2K(Z(36))+,5*TCINC 

DTMINsZ(37)/CKFT 

DRMINsZ (38) 


D-40 



Astronuclear 

Laboratory 


W1INC«Z(39) 

K0INOZ(«0) 

RlMAXaQ. 

IF (KIINC.GT.O.) R1HAXsZ(41)+0.5*RIINC 

R9MAXS0, 

IF (ROINC.GT.O . ) R9MAXsZ(42)+0.5*R0INC 
PWRDENsZ (43) 

UINPC«Z (44) 

OGAM*Z (45) 

RLPERT*Z(46)/PCMULT 

RLMAX*Z(47)/PCMULT 

IF (RLMAX.NE.O.) NOPTIM«l 

ZlIBRHsZ(50)/ZLl 

ZLIBRC«Z(5l)/ZLl 

l K9*Z(52) 

CTEIC«Z(53)*1.0£-06 
CTE0C = Z(54)M .0E-06 
DEL7L*Z (55) 

HRINC«Z(56) 

HRSslOOO. 

IF (HRINC.LE.O.) GO TO 5 
HRSsQ. 

irite=o 

5 HRHAXs Z (57) 

HFLF®Z(58)/7,907E-05 

TREJ*Fk2K(Z(59)) 

RADPCT aZ (60) 

WTCON«Z(fel) 

OGRFsZ (62) 

DGTRFapK2K(Z(63)) 

DGLRFa.01*DGRF*Z(64) 

RADKaZ (65) /CKFT 
C IDP«Z (66) 

RADCa50l.E-l2*Z(67)*Z(68)*Z(69) 

DO 6 I a 1 » 30 

6 DP(I)b2(70+I) 

IF (KA3E1.E0.1.0R. (IRITE,NE.3.aND.IRITE.NE.4) ) 

1 CALL PHONYCDP, IPUN#HRS) 

KASEiaKASEl+I 
RLPC ( 1) aRLQAD 
YNCaZNC 

IF (NGT1.GE.3) YNCaPCHULT 
DO 97 lsl#NC 
97 RLPC(I)»RLOAO/YNC 
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ZCIRsYnOZLTE 

NGT2«2 

IF (RLQAD.EG.O.) NGT2*3 
IF (NGT1.LE.U0R.NGT1, EG. 3) NGT2*1 
IF (VReQ.GT.O. .AN0.NGT1 .GE.4) nGT2»« 

IF (NGTl.GE.fe) N g T 2 * 5 
Ns 1 

IF (IRITE.EQ.1.0R.IRITE.GT.2) QO TO 11 
N s 0 

C * WRITE INPUT NZ AND Z ARRAYS IF SPECIFIED 
WRITE (6,100) (CMT(J) , J*l,25) 

100 FORMAT (1HI,«X,2A4, IX# 18A4,2H ( , F 2 .0, 2 ( lH/ , F2 , 0 ) , 2H) ,2A«) 

WRITE (6,101) 

101 FORMAT (1H0/50X, 10HINPUT DATA/) 

DO 103 I«l,10 

103 WRITE (6, 10a) (J,NZ(J) ,J=I,50,10) 

104 FORMAT (2X f 5(3HNZ(,I2,2H)»,I4,llX)) 

WRITE (6,112) 

DO 105 1*1,20 

105 WRITE (6,106) (J,Z(J) ,J«I, 100,20) 

106 FORMAT (2X,4(3H Z ( , 12 , 2H ) * , G 1 2 . 4, 3X ) , 2hZ ( , I 3 , 2H ) « , G 1 2 , U ) 

CHECK TO SEE HOW TEMPERATURES ARE TO BE READ IN 

11 IF (NGT1.GE.3) GO TO 5080 

WRITE (6,100) (CMT(J) , Jsl,25) 

IF (INT.NE.1) GO TO 5070 
WRITE (6,5069) 

5069 FORMAT (1H0,25X,61H***** SAME AXIAL TEMPERATURES USED AS IN PREVI 
10US CASE *****) 

K F T E M s i 
GO TO 123 

5070 TH A V * 0 , 

TCAVsO, 

IF (INT.NE.O) GU TO 5072 

READ (5,5071) (TEMP(l,I),TEMP(q,I),Isl,Nc) 

5071 FORMAT (12F6.2) 

IF (1RITE.LT.3) WRITE (6,5075) TEMKF 
5075 FORMAT (IH0,22X,62HH0T ANO COLD CLAD TEMPERATURES OR EACH AXIAL SE 
ICTION IN DEGS, , A 1 ) 

GO TO 123 

CALCULATION OF EQUI-SPACEO (STANDARD) THERMOCOUPLE LOCATIONS 

5072 IF (INTTC.NE.O) GO TO 5083 
DTAVsZcIR/(FLOaT (INT)-U) 

XTHERM ( 1 ) *0 , 

00 5082 1*2, INT 
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5082 XTHERM(I)aXTHe«M(I-l)+DTAV 
GO TO 5 084 

5083 READ (5,5071) ( X THERM ( I ) , I a 1 , InT ) 

5084 READ (5,5071) (THERMH(I) , 1*1 » InT) 

READ (5,5071) ( THERMC ( I ) , I ■ 1 , InT) 

IF (ZUP ,GT, 0.) GO TO 67 

WRITE (6,66) 

66 FORMAT (/////, 7X f 96HTEMPERATURES CAN NOT BE INTERPOLATED ON A RUN 
lOPTlMUING AN AXIAL LENGTH, CALCULATIONS SUPRES8ED) 

GO TO 93 

CHECK FOR ANO ELIMINATE ZERO THERMOCOUPLE INPUT TEMPERATURES 

67 J = 0 
K *0 

DO 65 1*1, INT 

IF (THERMH(I) ,LE,0,) GO TO 63 
J*J + 1 

VPC(J)*THERMH(I) 

VQC(J) aXTHERM(I) 

63 IF (THERMC(I) .LE.O.) GO TO 65 
K*KM 

PE (K ) * THERMO ( I ) 

RPC(K)*XTHERM(I) 

65 CONTINUE 

UTAVs,5*ZLTE 
00 81 1*1, NC 

TEMPC1, I)*SI(VOC, VPC,DTAV, J, IND,NXTRP,NXTRP) 

TEMP (9, I)*SI (RPC, PE#DTAV,K, IND,NXTRP,NXTRP) 

81 dtav*dtav+zlte 

WRITfc (6,5074) TEMKF 

5074 FORMAT (IH0,27X,51HINPUT THERMOCOUPLE LOCATIONS AND READINGS IN DE 
1GS, , A 1 ) 

NKN*(lNT+9)/10 
DO 73 K* 1 , NKN 
KNK *1 0*K 
KNNK*KnK-9 

IF (KNK, GT, INT) KNK*INT 
WRITE (6,108) (I, I«KNNK,KNK) 

WRITE (6,112) (XTHERM(I),I*KNNK,KNK) 

112 FORMAT (1X,F8,3,9F1 1 ,3) 

WRITE (6,112) (THERMH(I) , I»KNNK,KNK) 

73 WRITE (6 » 112) (THERMC (I ) , I aKNNK , KNK ) 

IF (IRlTE.GE.3) GO TO 122 
WRITE (6,5076) TEMKF 

5076 FORMAT ( IH0/2 1 X , 67HI NTERPOL ATEO MID-POINT TEMPERATURES OR EACH AXI 
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1AL SECTION IN DECS. ,A1) 

123 NKNs(Nc+ 9)/10 
DO 111 K«1,NKN 
KNK«lO*K 
KNNKaKNK»9 

IF (KNK.GT.NC) KNKaNC 
WRITE (6,108) (I, I*KNNK,KNK) 

108 FORMAT (1HQ.3X, 10(13, 8X) ) 

WRITE (6,112) (TEMP(1 , I) , I»KNNk,KNK) 

111 WRITE (6,112) (TEMP(9,I),I«KNNK,KNK) 

1 22 IF (INT.eO.l) GO TO 5073 
DO 5077 I * 1 , NC 
TEMR(l,l)sFK2K(TEMP(l,I)) 

TEMP (9, 1)nFK2K(lEMP(9,D) 

THAVsTHAV + TEMPd, I) 

507 1 TCAV»TCAV*TEMP(9 # I) 

THAVaTHAV/2NC 

TCAV*tTCAV/ZNC 

DTAVsThAV-TCAV 

TM0DAVsTCAV+0,5*DTAV 

THAVF«CKF2(2)*THAV-CKF1(2) 

TCAVF=CKF2(2)*TCAV-CKF1 (2) 

TM0DFacKF2(2)*TM0DAV-CKFl (2) 

DTAVF«CKF2(2)*0TaV 

5073 WRITE (6,5078) Tha V , TH A VF , TC A V , TC A VF , DT A V , OT A VF , T MOD A V , TMODF 
5078 FORMAT ( 1H0/22X, 30HAVERAGE HOT CLAD TEMPERATURE a,F9.3, 9H DEG.K 
1( ,F9.3,7H DEG,F) //21X, 31HAVERAGE COLD CLAD TEMPERATURf *,F9.3, 

2 9H DEG.K (,F9.3,7H DEG . F ) // 1 BX , 34HA VE&AGE RADIAL TEMPERATURE DR 
3UP s,F9,3,9H DEG.K ( , F9 , 3, 7H DEG.F) //24X , 28HA VERAGE MODULE TEMPER 
4ATURE a,F9.3,9H DEG.K (,F9.3,7H DEG.F)) 

GO TO 119 

5080 IF (INT.EO.l) GO TO 119 

READ (5» 507 t ) (TEMP(l,I),TEMP(9,I),Ial,NC) 

IF (N.FO.O) WRITE (6,62) TEMKF , ( TEMP ( 1 , I ) , TEMP (9, 1 ) , I a 1 , NC ) 

62 FORMAT (1H0/36H INPUT TEMPERATURE PAIRS (IN DEGS. ,A1,1H)/ 

1 (6(2X,F'7.2, 1H-,F7.2, 1H, ) ) ) 

IF (KFTEM.EO.l) GO TO 119 
DO 5081 1*1 , NC 
TEMP(1,I)«FK2K(TEMP(1,I)) 

5081 TEMP(9,I)sFK2K(TEMP(9,I)) 

119 IF (RLOAD.GE.O.) GO TO 118 

READ (5,5071) (RLPC ( I) , I ■ 1 , NC) 

IF (N.fcQ.O) WRITE (6,117) (KLPC(l)#Ial,Nc) 

117 FORMAT (74H INPUT LOAD RESISTANCE FOR EACH TEMPERATURE PAIR (IN MI 
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lLLIOHMS PER MODULE) / (2X, 1 0 CGI 0 . 2 » l H, ) ) ) 

DO lib I ■ 1 # NC 

lib RLPC(I)a.001*RLPC(I)/PCMULT 
lta DO 125 IalfNC 

DTMOO(l)s*TEMP(l» I) -TEMP (9# I) 

IF (DTMOD(I) .Nt.O,) GO TO 13b 
TEMP (1,1) ■TEMP(1, I) ♦.001 
DTMOD(I) * .001 

12b IF((TEMP(l#I).GT.0..AND.TEMP(9,I).GT.0.).0R,NGTl.GE.3) GO TO 12b 

WRITE (b, 121 ) 

121 FORMAT (1H0,/17X,80HNEGATIVE OR ZERO TEMPERATURE ENCOUNTERED - CAl 
1CULATI0N SUPRESSED FOR THIS CASE ) 

GO TO 93 
12b IT(I)*0 
127 CONTINUE 

C ” ESTABLISH MODULE RADIAL DIMENSIONS FROM INPUT DATA 
3061 IF (Z(2).LE.O.) DR ( 1 ) =C I 02*H AO ( 1 ) +C IDP 
RAD(2)»RAD(l)+DR(l) 

RAD(3)*RAD(2)+DR(2) 

RAD(a) B RAD(3)+DR(3) 

RAD(5)*RAD(9)+DR(9) 

IF UZ9.E0.1) GO TO 3068 
DRC5J8Z9SAV 

3059 RAD(b)sRAD(5)+DR(5) 

RAD(7) sRAD(6)+DR(b) 

RAD(8)«RAD(7)+DR(7) 

IF (DR8.LE.0.) GO TO 3067 
DR(B)sDR8 

RAD(9)sRAD(8)+DR(8) 

GO TO 3066 

3067 IF (C0D2.EQ.1.) GO TO 1996 
KAD(9)aRA0(8)/(1.0-C0D2) 

DR(8)»RA0(9)-RAD(8) 

GO TO 3066 

3068 RAD(9)aZ9SAV 
3070 DR(8 )»dR8 

IF (DR8.LE.0.) DR(8)*C0D2*RAD(9) 

RAD(8)»RAD(9)-DR(8) 

RAD(7)aRAD(8)-DR(7) 

RAD(6)*RAD(7)-D«(6) 

DR(5)3RAD(6)-RAD(5) 

3066 Z0D»2 ,0*RAD (9) 

ZID«2.0*RAD(1) 

IF (RRREF.LE.O.) RRREFsZOD/ZlD 
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I H NOK* J 

IF CENOZK.GT.O. ) IENDK*2 
0G£NL*5l ,48*PWR0EN*RA0(1) *RAD(1 ) 

UGENsQGENL*ZLTE 

IF (IOlJMZ.NE.l) GGEN = QGENL*2,*(ZLN4-ZLI) 

IF (QInPC.GT.O,) QGENbQINPC 
IF (IRlTE.GT.3) GO TO 4004 

C - WRITE RADIAL AND AXIAL DIMENSIONS IF SPECIFIED (NZ(5) ,FQ, 0,1,2) 
WRITE (6,100) (CMTtJ)# J«l,25) 

HD1 Cl)aSSlN(NCLOH) 

HDl (3 )sFETm(ICRH) 

HD1 (4)«CASCAD(l3N,2) 

HD(1,4)sCASCA0(I3P,2) 

HDl (5) «CASCAD(I2N,2) 

HD(1,5)*CASCAD(I2P,2) 

HDl (6)*FETM(ICRC) 

HDl («)*SSIN(NCLUC) 

WRITE (6,2017) 

2017 FORMAT (1H0//, 10X, 17HRADIAL D IMENSIONS, 1 3X , 1 8HR ADI AL THICKNESSES, 

1 7 X , 9HC0MP0NENT,/) 

1*1 

IF (NOPTIM.NE.O.AND.DRINCR.GT.O.) 1*2 
DO 2019 J* 1 , 8 

2019 WRITE (6,2018) J , R AD ( J ) , HD 3 ( J , I ) , J , DR ( J ) , HD3 ( J + 8 , I ) , HD 1 ( J ) , MD2 ( J ) , 
1 (HD(K, J),K«i,6) 

2018 FORMAT (10X,4HRAD(, II, 3H) *,F10.6,5M IN, , A 1 /40X , 3HDR C , 1 1 , 3H) ■, 

1 F10,6 f 4H IN.,Al,3X,A3,Al,6A4) 

J *9 

WRITE (6,2018) J, RAD(J) 

J 1 • 1 

IF CNOPTIM.NE.O.AND.VREQ.GT.O, .AND.(ZLREQ.GT.0..0R.PEREQ.GT.0.) ) 

1 J l *2 

J * ( J 1 ♦ IDUMZ+ 1 ) /2 

WRITE (6,2021) ZLN,ASTAR1 (Jl) ,ZLP,ASTAR1 (J) , ZL I , ZLTE , ASTAR1 ( J) 

2021 FORMAT ( 1 HO/ /46X , 1 8HAX I AL DIMENSIONS , / / / 1 OX , 8HZLN *,F9.5, 

1 23H IN, (N-TYPE WASHER) , A2//10X, 8HZLP »,F9.5,23H IN. IP- 
2TYPE WASHER) A2//10X,8HZLI «,F9.5,24H IN, (MICA THICKNESS)// 

3 10X,8hZLTE * , F 9 , 5 , 32H IN, (TOTAL COUPLE THICKNESS) , A2) 

IF (J.EO.l) WRITE (6,2023) ZClR,YNC 

2023 FORMAT (1H0, 9X,8HZCIR »,F9,5»30H IN, (TOTAL CIRCUIT LENGTH,, 
1 F4,0,9H COUPLES)) 

IF (J.EQ.2.0R.I.EQ.2) WRITE (6,2022) 

2022 FORMAT ( 1 HO, 50X , 5 1 H* THIS DIMENSION OPTIMIZED OR CALCULATED INTERN 
1 ALL Y ) 
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4004 DR4sQR{4) 

0RTE»RAD(b)-RA0(4) 

IF (ORTE.LT.DRMlN.AND. (NOPTIM.GT.O.OR.NCLAD.GT.O) ) 60 TO 199b 
3016 IF (RAO(l) .LE.O. .OR. ZLI, LE.O.) GO TO 199b 
DO 31 J* 1,8 
J1*J+1 

IF (RAD(Jl).LE. 0.) 60 TO 199b 
LNRAD(J)*AL0G(RAD(J1) /RAD(J)) 

31 IF (J/2.NE.2,AND,LNRA0(J) .LE.O.) GO TO 1991 
I SEC T a l 
GO TO 1999 

1991 CONTINUE 

IF (NOPTIM.EO.O.OR.ROINC.EO.O.) GO TO 1996 
ISEC T *3 
GO TO 1999 
1996 WRITE (6,1992) 

1992 FORMAT (1H0,21X,68HERRUR IN DIMENSIONS SPECIFIED - CALCULATIONS SU 
1 PRESSED FOR THIS CASE) 

GO TO 93 

1999 ZKEN0«ENDZK*AL0G(RRREF)/AL0G(Z0D/ZID) 

IF (NOPTIM.EQ.O) GO TO 4002 
CALL OPTIMUSECT, J) 

GO TO (4002,4002,4002,3059,4002,3070,3061), J 
CALCULATE ALL RECURRING PRODUCTS 
4002 ALR64 s LNRAD (5) +LNRAD £4) 

CGAMa,5/ALR64-l,/(EXP(2.*ALR64)*l.) 

CONST l*. 062659* ALR64 
CONRN«cONSTl/ZLN 

RC0N4 ■ REC0NH/(40,5366*ZLN*RAD(4) ) 

RCONb s REC0NC/£40.5366*ZLN*RAD(6)) 

CHECK FOR THERMAL CONTACT COEFFICIENTS EQUAL ZERO 
CONST 3a0 , 

CONST 4o0 , 

IF (HH.NE.O.) C0NST3a t 3937*(1.0/RAD(2)+1.0/RAD(3))/(LNRAD(2)*HH) 

IF (HC.NE.O.) CONST4a.3937*(1.0/RAD(7)-*-l,0/RAD(8))/(LNRAO(7)*HC) 
CHECK TO SEE IF P»LEG AXIAL LENGTH IS TO BE OPTIMIZED 
IF (IDUMZ.GT.l) GO TO 78 

CALCULATE RECURRING PRODUCTS IF P-LEG IS NOT To BE OPTIMIZED 
ZLPNaZLP+ZLN 

zlpni*zlpn*zli 

zlte*zlpni+zli 

zlnlp«zln/zlp 

C0NST5al5,9593*ZLP 

ZLIP»2.0*ZLI/ZLP 
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C0NST2b15.959.WlTE 

HINT4srC0N4*ZLPN//lP 

WInT6sRCUN6*ZLPN/ZLP 

CONWPscONST 1/ZLP 

IF (ICPMPL .NE. 0) GO TO 4005 

ClFEM*l5.9593*DR(3)*(RAD(3)+RAn(4))/ZLTE 

ClFEC*l5.9593*DR(6j*(RA0(S)*RAn(7))/ZLT£ 

GO TO 78 

CALCULATE RADIUS OF CIRCULAR CONDUCTOR RING ELECTRICAL STREAMLINE 
4005 HCHOTsaMINI ( DR C 3 ) , ZLN , ll P } 

KCOLDsaMINI COR (6) ,ZLN,ZLP) 

C1FEH=5»0S*RAD(4)*ALUG(3,14159*RCH0T/ZlI*1,0) 
ClFEC*5,Oft*RAO(6)*ALOG(3.14159*RcnLD/ZLl+l ,0) 

78 IN()I = InOI + 1 
C - ENTER SUBROUTINE COUPLE 
CALL COUPLE 

IF (NO n U MP.EU.5J GO TO 93 
C - ENTER SUBROUTINE LIFE, IF SPECIFIED 
IF (HRINC.LE.O) GO TO 3052 
CALL LIFE(DAPERI) 

IF (J0PTIM.LT. 01 GO TO 93 
GO TO 78 

3052 IF (NOpriM.EQ.Ol GO TO 3051 

C - ENTER SUBROUTINE OPTIM FUR THOT, TCOLD, PCT2N PERTURBATION 
CALL UpT I M (2 , J ) 

OTMOOClJsTEMPd, II -TEMP (9, 1) 

GO TO (93, 4002,78,3059, 78 , 3070, 3061 ) , J 
3051 IF (NPIJMP.LE.O) GO TO 3053 
CALL PUMP(J) 

GO TO (4002 , 30b 1 , 93) , J 

3053 IF (NOoUMP.EQ.3J GO TO 93 
IF (NGT1.GE.3J GO TO 79 

CALCULATE MODULE VOLTAGE AND RESISTANCE FOR CURRENT CONVERGENCE CHECK 
UC V*0. 

QHM*0 . 

UMOD=0. 

00 35 1*1, NC 
UCVsOCV+VOC (I J 
OHM*OHM + RPC (I J 
3b QMOD*QmOO»QT Cl) 

J * I NO I + 1 

IF (NGT2-2J 79,37,22 
22 RLUADsOHM 

CHECK FOR CONVERGENCE ON MODULE CURRENT 
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37 CURMOD(J) « OCV/tOHM+RlOAD) 

IF (ABS(1.0-CUR«0D(INDI) /CURMODCJ) ) .LE , 1' OLCUR . AND , I NO I . GT .2) GOTU79 

IF (INDI.EU.2I T0LTEM«0.1*TnLTEH 

IF (INDI .LT, MAXIND) GO TO 78 

WRITE (6. 43) MAXIND, ( J , CURMQD ( J) , J * 1 , IND I ) 

43 FORMAT (///,27X»37H***CURRENT CLOSURE NUT OBTAINED AFTER, IS, 14h 
HTERAT l0NS***,//,22X,b5HM0DULE CURRENT CALCULATED ON EACH PASS I HR 
2UUGH SUBROUTINE COUPLE , // / , 42X , 9HI TER AT ION, 7X , 7HCURRENT , ✓ , 
3(IH0,44X,I3, IH.,7X,1PE12.5,/)) 

74 CALL RITE 

IF CNGT1.EQ.5.AND.NGT2.E0.3) GO TO 2060 
IF (NGT1.NE.0) GO TO 2059 
I NO I * 1 
NGTU2 

toltem«io.*toltem 
IF(RL 0A0.GT t 0.) GO TO 2053 
CURMUD(2)«0.50*0CV/0HM 
NGT2*3 
GO TO 78 

2053 NGT2b 2 
RLPC(l)aRLOAO/ZNC 

00 2054 I ■ 1 , NC 

2054 RLPC (i)aRLPC(l) 

CURM0D(2)e0CV/(0HM+RL0A0) 

GO TO 78 

2059 IF (NGT1.NE.5) GO TO 2061 

2060 CURM0D(2)«0. 

1 NO I * l 
NGT2«1 
NGTU3 

dttesTemp( 4, n - temp ( 6, n 
GO TO 78 

2061 IF (NGTI.NE.7) GO TO 93 
NGT 1 *4 

IRITE*4 
INTal 
GO TO 32 
END 
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SUBROUTINE RITE 

DIMENSION TAVG(300) ,ETAC(300) ,0TMF(300) , TF(9, 300) , T I CIJP (300) , 

1 nTF(9,300) ,CPMD( 12,2) r ALPCUP( 300) , ETA (300) 
COMMON/TERIT/DTAy,OHM,OCV,RLOAD>QMOD,TCAV 

COMMON /M0CP0P/RA0(9) #0R(8) , TEMP (9 , 300 ) » VOC (300) ,RPC(300) » PE (300) , 

1 QT(300) ,NOPTIM, VPC(300) ,RLPC(300) ,DTM0D(300) ,CUR(300) ,LNRAD(9) , 

2 ZlP,ZlN,ZLI, ZLNLP,ZLPNI,ZLTE,NG71,NGT2, ITPERT,PCMUI.T,NODUMP, IVD, 

3 ogen,ogenl,qte*qout,zid,zoo,zkend,zk9,2kr,rppc»rpn, toi,tem,trej> 

4 vkeu,pereq, zlweq, wtcon, zlpn, idumz,rn, dgrf,dgtrf,dglrf,radk, jdump, 

5 DP (30) , NC 

COMMUN /M0DCPL/CURMQ0(15) , IT (300) >OTTE»CONRN,CONRP,ZLIP» znc, 

1 ALR64, MAXTEM,QGAM,CGAM, INDI,RINT4,RINT6,RC0N4,RCDN6,C1FEH,C1FEC, 

2 CONST 1 , CONST 2, CONST 3, CONS T4,C0NST5, NR50pT,SEGTE, ICRMDuR ADC 
COMMON /MOOOPT/IPOWT, 0R4 . DR INcR * I« I TE , DWT E » JOPT IM , I Z9 , IENQK . IPUn , 

1 RI lNC,ROlNCr R1MAX,R9MAX,DTMIN,QREQ, TH I Nc , TCINC , TCMAX , THM AX 
it NPFRr,ZNPERT,CKFP,CKFT,CT£0C#cTElC»DELTL,RLPERT,RLMAX,ZNMlN,Z9SAV 
COMMON /TITUE/FETMC3) ,CASCAD(8,2)»SSIN(3) ,CMT(25),C0N0(3,5) ,TEMKF, 
1 CRMD(2,2),I3N,I2N,I3P,12P 
COMMON /CPLRIT/DT(8,300),QPTZLP(300) ,0PTCPL(300) 

DATA CPM0/4H C IR , 4HCU I T , 4H C A , 4HRN0T , 4H CUR * 4HRENT , 4H EFF,4H (AM, 

1 4HPS) i 4H PCT,1H.,4H. ,4H MO, 4HDULE r 4H A VE , 4HRAGE t 4H R-L, 

2 «HOAD ,4HTEMP,4H(M-0,4HHMS),4HDEG.,4HMDHM, 1HS/ 

CPM0(1,2)=CRMD(1, 1) 

CPMD(2,2)*CRM 0(2,1) 

WRITE (6,100) (CMT(J) , Jsi ,25) 

100 FORMAT (1H1, 4X,2A4, IX, 18A4, 2H ( , F2 , 0 , 2 ( 1 H/ , F2 , 0 ) , 2H ) ,2A4) 

IF (ZKFND.LE.O.) GO TO 11 
CPMD(1,2)*CRMD(1,2) 

CPM0(2,2)«CRMD(2,2) 

11 IF (NGT1 .LT.3,aN0.IR1TE,GT.2) GO TO 200* 

DO 2001 1*1, NC 

P* TEMP (9, I)+0,5*DTMQ0(I) 

TAVG(I)*0K2FK(P) 

ETACCI)«100,*DTM0D(I)/TEMP(1,I) 

DTMF(I)*CKFT*0TM0D(I) 

TF(9,I)*DK2FK(TEMP(9,I)) 

ti(:up(i)*dtmf (i)/ot(i) 

ALPCUP(I)*VOC(I)/DTMF(I,) 

ETA(I)*P£(I)/0T(I) 

DO 2001 J * 1 , 8 

DTP (J, I)*CKFT*DT (J,I) 

2001 TF(J,I)»DK2FK(TEMP(J,D) 

200b IF (IRITE.GT.2) GO TO 2007 
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1994 WRITE (6,1998) ( COND (.1 , NGT2 ) , Jsl , 3) , T EMKF , ( J , J« 1 , 9) , 

1 (I, (TF(J»I),Jsl,9)»I=l»NC) 

1998 FORMAT (1H0, 1 OX , 58HR AD I AL TEMPERATURE PROFILE FOR EACH COUPLE CALC 
1ULATED FOR ,3A4,19H CONDITIONS (PEGS, ,A1,1H)/8H0 CQUPlE,6X, 

2 9(1HT, I1,9X)/(3H0 , 1 3 , 1 H . , F 1 1 .2 , 8F 11 . 2) ) 

IF (NC.GT.7) WRITF (6,100) (CMT ( J) , J« 1 , 25) 

WRITE (6,2009) TEMKF , (I , (OTF ( J , I ) , J« 1 , 8 ) , I s 1 , NC ) 

2009 FORMAT (1H0,21X,64HRADIAL TEMPERATURE DROP ACROSS COMPONENTS FOR E 
IACH COUPLE (DEG. ,A1,1H)//8H COUPLE , 1 OX , 28HINNER BORON 

2 INNER, 28X,26H0UTER BORON OUTER/ 1 9X , 81 HCL AO NITRIDE 

3 CONDUCTOR T/El T/E2 CONDUCTOR NITRIDE CLAO/ 

4(1H0,I5,2H, ,F16.3,7F11,3) ) 

WRITE (6,100) (CMT (J) » J* 1 , 25) 

2007 WRITE (6,2002) (COND ( J , NQT 2) , J t 1 , 3) 

2002 FORMAT (1H0, 12X , 61 HCALCULATEO PARAMETERS FOR INDIVIDUAL COUPLES OP 
1ERATING UNDER ,3a«,11H CONDITIONS) 

C - WRITE INDIVIDUAL COUPLE PARAMETERS 

WRITE (6,13) (CPMD(J,l),Jal, 7), TEMKF, CPMD(8,1),CPMD(9,1), TEMKF, 

1 TEMKF, CPM0(10,1),CPMD(11,1) 

13 FDRMAT (1H0,12X,18H*****V0L'TAGeS****M2X,50HINTERNAL POWER 

1 HEAT THERMAL EFFECTIVE , 2 a4, 1 X , 2A4/33H TH / TC VOC 

2 V-LOAD ,2A4,63H RESISTANCE OUTPUT INPUT IMPEDANCE SEEbE 

3CK EFF. , A4, 1H./8H (DEG. , A1 , 24H) (VOLTS) (VOLTS) , 

4 2A4,35H (M-OHMS) (WATTS) (WATTS) (D-, A 1 , 1 2H/K W) (MV/D*, 

5 A 1 , 1 4H ) ( P C T . ) ( , A4 , A 1 , 1H) ) 

WRITE (6,14) (TF(1,I),TF(9,I) , VOC(I), VPC(I),CUR(I),RPC(I),PE(I)r 
1 QT C I ) ,TICUP(I),ALPCUP(I),ETA(I),ETACCI) , 1=1, NC) 

14 FORMAT (1H0,F5.0, 1 H/ , F4 , 0 , 2F 10, 5, F 1 0 , 3, 3PF 1 0. 4 , OPF 1 0. 4 , F 1 0, 2 , 

1 3PF10.2,3PF10.4,2PF10,3,OPF9.2) 

C - WRITE CALCULATED AXIAL LENGTHS IF ZLP HAS BEEN OPTIMIZED 
2062 IF (IDUMZ.EQ.1) GO TO 2061 

IF' (NC.GT.7) WRITE (6,100) (CMT ( J) , J*1 , *5) 

WRITE (6,2052) ( I , ZLN, OPTZLP C I) , ZL I , OPTCPL U ) , OPTCPL ( 1 + 1 00) , I a 1 , NC ) 
2052 FORMAT (1H0// 14X, 6HC0UPLE , 1 8X , 52HAX I AL THICKNESSES CALCULATED I 

IN OPTIMIZATION R0UTINE/14X,6HNUMBER, 1 1X,5HN-LEG, 10X , 5HP-LEG, 7X , 12H 
2MICA W ASHERS, 6X, 6HC0UPLE, 8X, 7HCIRCUIT/ (1 HO, 13X, 14, 4H) ,5F15.5) ) 

2061 IF (NGTt.LT. 3) GO TO 2005 

IF (NC.GT.7. OR. IDUMZ.EO, 2) WRITE(6,100) ( CMT ( J ) , J* 1 , 25) 

WRITE (6,2042) PCMULT, (COND ( J , NGT2) , J= 1 , 3) 

2042 FORMAT (// F7 • 0 , 32H COUPLE MODULES OPERATING UNDER ,3A4,S5H CUND 

1ITI0NS WITH UNIFORM HOT AND COLD CLAD TEMPERATURES) 

RLD«0, 

WRITE (6,13) (CPMD(J, 2), J*l, 7), TEMKF, CPM0(8, 2), CPMD(9, 2), TEMKF, 

1 TEMKF, CPMDCIO, 2), TEMKF 
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DO 57 J a 1 , NC 

OCV « PCMULT*VOC(I) 

V LU = PCHULT*VPC(n 

IF (NGT2.NE.1) RLDslOOO,*PCMULT*WLFC(n 
OHMsPCmULT *RPC (I) 

P a PCMUIT*PE(I) 

QM0Q«PcHULT*QT'Cl)+ZKtND*nrHC10(l) 

sebspcmult*alpcup(i) 

f:taosP/omod 

TIUaDTMF(I)/UMOO 

57 WRITE (6,14) TFC1,I)|TF(9,I),OCV,VL0,RLO,OHM,P,QMOD,TIO,SEB,ETAU, 

1 TAVG(I) 

GO TO 2221 
2005 WRITE (6,2059) 

2059 FORMAT (1H0///43X,24H0VFRALL MODULE OPERATION) 

WRITE (6, 13) (CPM[)(J,2) , Ja 1 , 2) , ( CPMD ( J , 2) , J*5 , 6) , (CPMDfJ, l) , J»5,6) , 
lCPMOtlg, 1) ,TEMKF, (CPMD(J,1) , J«6,9),TEMKF,TEMKF, (CPMDCJ,2) , J«ll ,12) 
C - CALCULATE MODULE PARAMETERS 
rAVG(2)aOK2FK(TCAV) 

TAVG( l)aTAVG(2)+CKFT*DTAV 
VL D sOC V-CUHMOD ( 1N01 ) *OHM 
RLDelO0O.*VLD/CURM0DUNDI) 

PaVLO^CURMODCINDI) 

QMUDaQMOO+ZKEND*DTAV 

SEB*aCV/(OTAV*CKFT) 

ETAU*P/OMOO 

TIOsDTaV*CKFT/OMOD 

WRITE (6,14) TAVG(l) , T A VG ( 2) , OC V , VLD , CURMOD (I ND I ) , OHM , P, QMOD , T I 0 , 

1 SEB , E T AG , RLD 
WRITE (6,2011) I NO I 

2011 FORMAT (1HO,////,30X, 12,37N ITERATIONS TAKEN FOR CURRENT CLOSURE) 
2221 WRITE (6,2016) ( I T ( I J • I« 1 . NC) 

2016 FORMAT ( 1 HO// 38HONUM0ER OF ITERATIONS FOR EACH COUPLE-/ (5X, 25 14) / 

1 1H0#22(5H*****) ) 

RETURN 

END 
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SUBROUTINE PHONY (DP, IPUN,HRS) 

DIMENSION FM3N(15) ,FM2N(lS) ,FM3P(15) ,FM2P(15) ,DP(30) 

COMMUN /TITLE/FETM(3),CASCAD(8,2),SSIN(3),CMT (25) , CONO ( 3 , 5) , T EMKF , 
1 CRMD(2,2) , I3N, I2N, I3P, I2P 

COMMON /DATIN/TTT(15) ,ZKMIC1 (15) # ZKMlc;2(15) , ZKSS ( 1 5) , ZKBN ( 1 6 ) , 
l ZKIN(15) ,Z2NA(15) ,Z2NR(15) ,Z2NK (15) ,Z2PA(15) ,Z2PR(15) ,Z2PK(15) , 
2Z2NRK(l5),Z2PRK(15) ,Z3NA(15) , Z3NR ( 15) , Z3NK ( 1 5) , Z3NRK ( 1 5) , Z3P A ( 1 b) , 
3 Z3PR(l5),Z3PK(15),Z3PRK(15),ZKCRH(15),ZKCRC(15),RCRH(l5),RCRC(l5) 
4, RHP(15) rRPRHO 

COMMON /F0DAT/ALPHA(15,8) , RHO ( 1 5 , 8 ) , ZK ON ( 1 6 , 8 ) * ZKSS I N ( 1 6 , 3 ) , 

1 RFETM(15,4) ,ZKFETM(lb,3) ,ZKMIcA(lb) 

COMMON /F0wT/DEN(8) ,OENI, DEN3N,OEN2N, 0EN3P,DEN2P,DCRH,DCRC 
COMMQN/TEFO/NCLOC* NCLOH# ICRH, IcRC , I P I N , Z l I BRH , ZL IBRC , TP I N , RP I N 
DATA H0T/4HH0T / , C0LD/4HC0LD/ , NNAME/2HN-/ , PNAME/2HP-/ 

IF (NCLOH. NE,1. AND. NCLOH. NE. 3) NCLDH*2 
IF (NCLDC.NE.2,AND.NCLDC.NE.3) NCLDC«l 
IF (ICRH.NE.2.AND.ICRH.NE.1) IC«H*3 
IF (ICRC. NE, 2. AND. ICRC. NE. I) IC*C*3 
IF (I3N.LT. 10) GO TO 301 
I3N*I3 n-10 

READ (5,3001) (ALPHA(I,I3N),I*1,15), (RH0(I.I3N),I*1,15), 

1 ( Z K 0 N ( I , I3N) , 1*1, 16) ,CA3CA0(I3N, 1) , CASCAD ( I3N, 2) 

3001 FORMAT ( 8F9 . 2 /7F9 . 2/8F9 . 2/ 7 F9 . 2/8F9 . 2/8F 9 . 2 , 1 X , 2A4 ) 

IF (IPUN.GT.O) WRITE (7,401) (C ASC AD ( I 3N , I ) , I • 1 , 2 ) , ( AL PHA ( I , I 3N) , 

1 1*1,15), (CASCAQ(I3N,I) , 1=1,2) , (RHO (I, I3N) , I » 1 # 15) , 

2 (CASCAD(I3N, I), 1*1,2), (ZKON ( I , I3N) , I * 1 , 1 6 ) 

401 FORMAT (25HC SEEBECK COEFFICIENT OF ,2A«,23H MATERIAL (VOLTS/DEG. K 
1) ,5X, IhI, 10X,5(E9.3, 1H, )/5X, 1H2, 10X,5(E9.3, 1H, )/5X, 1H3, 10X,5(E9.3, 

2 1 H , ) / 1 7HC RESISTIVITY OF ,2A4,|9H MATERIAL (OHM*CM . ) /5X , 1 H 1 , 1 OX , 

3 5(E9,3, 1H,)/5X, 1H2, 10X,5(E9.3, 1H,)/5X, 1H3, !0X,5(E9,3, 1H, ) / 

4 26HC THERMAL CONDUCTIVITY OF ,2A4,27H MATERIAL (WATTS/CM. /DEG. K) / 

5 5X, 1H110X,5(E9.3, 1H, )/5X, 1H210X,5(E9.3, 1H, )/5X, 1H3,6(F9.3» 1H, )) 

301 IF (I2N.LT. 10) GO TO 302 
I2N3I2N-10 

READ (5,3001) ( ALPHA ( I , I2N) , I* 1 , 1 5) , (RHO ( I , I2N) , 1*1 , 1 5) , 

1 (ZKON(I, I2N),I« 1,1b), CASCAD (I2N, 1), CASCAD (I2N, 2) 

IF (IPUN.GT.O) WRITE (7,401) (c ASC AD ( I2N , I ) , 1 9 1 , 2 ) , ( ALPHA ( I , I2N) , 

1 I » 1 , 15) , (CASCAD(I2N,1), 1*1,2), (RHO ( I , I2N) ,1*1,15), 

2 (CASCAD(I2N, I) , 1*1,2) , (ZKON ( 1 , 1 2N) , I * 1 , 1 b ) 

302 IF (I3P.LT. 10) GO TO 303 
I3PSI3P-10 

READ (5,3001) ( ALPHA ( I , I3P ) , I* 1 , 1 5) , (RHO ( I , I3P ) , I * 1 , 1 5) , 

1 (ZK0N(I,I3P),I?1, lb) , CASCAD ( I 3P, 1) , CASCAD ( I3P, 2) 
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IF UPUN.GT.O) WRITE (7,401) (C ASCAD ( I 3P, I ) , I * 1 , 2) , ( ALPHA ( 1 , 1 3P ) , 

1 I = (CASCADCI3P, II # 1=1,2) # (RHO (I, I3P) #1 = 1# 15) , 

2 (CASCA0U3P, I) # 1=1,2) , (ZKON (I # I3P) ,1*1,16) 

303 IF (I2P.LT, 10) GO TO 304 
I 2P« I 2P* 1 0 

READ (5,3001) ( ALPHA ( I , I 2P ) , I * 1 # 1 5) , (RhO ( I , I2P) , I • 1 , 1 5) , 

1 ( ZKON ( I, I 2P) , 1*1, lb) ,CASCAD(I2P, 1 > , CASCAD ( I2P, 2) 

IF (IPljN.GT.O) WRITE (7,401) (c ASC AO ( 12P , I ) , £• 1 , 2) , ( ALPH A ( I # 1 2P ) , 

1 1=1,15), (CASCA0(I2P,I), 1=1,2) , C«MO ( I , I2P) , I = 1 , 1 5) , 

2 (CASCaD(I2P» I) , 1=1,2) , (ZKON (I , I2P) , 1*1 , 16) 

304 ZLIBHlsl .0-ZLIBRH 
ZLI8CU1.0-ZLIBRC 
PN3N * 1,0 

PN2N « 1,0 

IF (I3K1.EG.I3P) PN3NS-1.0 
IF (I2M.E0.I2P) PN2N*-1,0 
HRP=1,0E-S*HRS 
0/3NA*(1.+HRP*DP(I3n) )*Pn3n 
DZ 2NA* (1 ,+HRP*DP(I2N))*PN2N 
0Z3PA*1 ,0 + HRP*t)P (I3P) 

DZ2PA*1 ,0+HRP*0P(I2P) 

DZ3NR=1,0+HRP*DP(I3N+10) 

0Z2NR=1 .0 + HRP*j)P(I2N + 10) 

0Z3PR*l ,0+HRP*QP(I3P+10) 

QZ2PW*1 ,0+HRP*OP( I2P+10) 

HRP=,0l 

0Z3NK*1 ,U+HRP*0P(I3N+20) 

0Z2NK=i , 0+HRP*DP(I2Nt20) 

0Z3PK*1 ,0tHRP*DP(I3P+2O) 
t)Z2PK*l ,0 + HRP*uP{I2Pt20) 

1 = 16 

0EN(1 )=ZKSSIM(I,NCL0H) 

DEN(2)=ZK0N(I) 

l)EN(7)*ZKBN(I) 

0ENC8)*ZK8SIN(I,NCL0C) 

OENI«ZkNICA(I) 

OCRH»ZkFETM(I, ICRH) 

OCRC=ZkFETM(I, 1CRC) 

0EN3N*ZK0N(t, I3N) 

0EN2N*ZK0N(I, I2N) 

0EN3P*ZK(JN(I, I3P) 

DEN2P*ZK0NCI, I2P) 

DO 249 1=1,15 
RHP(I)*RFETM(I, IPIN) 
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ZKMICl (I) «ZLIBH1*ZKMICA(I)+ZLIBHH*2KFETM(I,2) 
ZKMlca(lJsZLIBCl*ZKMICA (I)+ZLIB«C*ZKFETM(I,2) 

RCRH(I)sRFETM(1, ICRH) 

RCRCCI)s«FETH(I» ICRC) 

ZKCRH(I WKFETMCI, ICRH) 

ZKCRC(I)«ZKFETH(I, ICRC) 

ZKIN(I)82kS3IN(I,NCLDH) 

ZKSscn *zks3in(i,ncldc) 

Z3NA(I)«DZ3NA*ALPHA(I, I3n) 

Z2NA(I)sDZ2NA*ALPHA(I, I2N) 

Z3PA(I)8DZ3PA*AUPHA(1, I3P) 

Z2PA(I)*DZ2PA*ALPHA(I, I2P) 

Z3NR(I)sOZ3NR*RHO(I»I3N) 

Z2NR(I)sDZ2NR*RH0(I, I2N) 

Z3PR(I)»DZ3PR*RH0(I,I3P) 

Z2PR(I)*DZ2PR*RHQ(I,I2P) 

Z3NK(I)»0Z3NK*ZKQN(I, I3N) 

Z2NK (I)aDZ2NK*ZKQN(I, I2N) 

Z3PK(I)»0Z3PK*ZKQN(I, I3PJ 
Z2PK (I)«DZ2PK*ZK0N(I, I2P) 

Z2NRK ( I ) s Z2NR ( I ) * Z2NK(I) 

Z2PRK ( i ) . Z2PR ( I ) * Z2PK ( I ) 

23NRK ( I ) a Z3NR ( 1 ) * Z3NK(I) 

249 Z3PRK(l)sZ3PR(I)*Z3PK(I) 

RPRHUaRPIN/SI(TTT,RHP,TPIN, 15,IN0,1,1) 
if cipun.lt. o) return 

WRITE (6,100) (CMT(J),J«1,25) 

100 FORMAT (1H1,4X,2A4, IX, 18A4,2H ( , F2 , 0 , 2 ( l H/ , F2 , 0) , 2H) ,2A4/ 

1 //38X,33HINPUT MATERIALS PROPERTIES TABLES) 

WRITE (6,281 ) HOT,NNAME, (C ASC AD ( I3N, I ) , Is 1 , 2) , HOT, PNAMf , (CASCAD 
1 (I 3P, I), 1*1, 2) 

281 FORMAT (1H0,/16X,2(A4,6H SIDE ,A2,13HLEG WASHER - ,2A4,9H MATERIAL 
l,3X)/l6X,2(42H**********************************«******* r 3X)) 

WRITE (6,282) 

282 FORMAT (8X,5HTEMP,,6X,2(39HSEEbECK RESISTIVITY CONDUCTIVITY 

1 , 6X ) ) 

WRITE (6,283) 

283 FORMAT (7X,6HDEG, X , 3X , 2 (45H (VOLT/DEG K) (OMM*CM) (WATT/C 

1M/K) ) ) 

WRITE (6,253) ( TTT ( I ) , Z3NA ( I ) , Z3NR ( I ) , ZiNK ( I ) , Z3PA ( I ) , Z3PR ( I ) , 
1Z3PKCI), 1*1,15) 

253 FORMAT ( 1H0, 6X , Fb . 1 , 8G 15 . 5 ) 

WRITE (6,100) (CMT ( J) , Ja 1 , 25) 

WRITE (6,281) COLD,NNAME, (C ASC AD ( I2N, I ) , I ■ 1 , 2) , COLD, PNaME , (CASCAD 
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1 (I2P,I),X«1,2) 

WRITE (6, 282) 

WRITE (6, 263) 

WRITE (6,253) ( T T T ( I ) , Z2NA ( I ) , Z2NR C I ) , Z2NK ( I ) , Z2PA ( I) , ?2PR ( I ) , 

1 Z2PK (I)»I*lfl5) 

WRITE (6,100) (CMT(J), J*l,25) 

WRITE (6,265) 

2 85 FORMAT UH0,/16X # a(12HFIG-0F-M(rRIT,3X) ,2(12HC0NDUCTIVITY ,3X) ,/8X, 

1 5MTEMP'.,2(7X.5HN-LES,10X,5HP-tE6,3X),3X # 27H8N INSULATOR INTER-C 
2UUPLE, /7X,6HOEG. k, 5X,2(9HM0T S I DE , 6X ) , 2 ( 9HC0L0 SIDE,6X),2aH SLEF 
3 V ES lNSULATORS/12X,a(6X,9H(i/DEG.K) ) , 1X,2(«X, V 1 H ( w AT T /CM/K ) ) ) 

00 250 1*1,15 

FM3N(I)cZ3NA (I) *Z3NA(I)/Z3NRK (I) 

FM2N(I)s*2na(I)*Z2NA(I)/Z2NRK(I) 

FM3P(I)sZ3pa(I)*Z3PA(T)/Z3PRK(i) 

250 FM2P(I)*Z2PA(I)*Z2PA(I)/Z2PRKCI) 

WRITE (6,253) ( 1 TT (I) ,FM3N(I) ,FM3P(I) ,FM2N(I) ,FM2P(I) ,ZKBN(I) , 

1 ZKMIC1 (I), 1=1, 15) 

WRITE (6,100) (CMT(J) , J=1 ,25) 

WRITE (6,286) HOT, COLD, SSIN(NClOH) ,SSIN(NCLOC) ,FETM (1CRH) , 

1 FETM(ICRH) ,FETM(ICRC) ,F£TH(1CRC) 

286 FORMAT (lH0,/17X f 25HINNEW CLAD OUTER Cl AD , 5X , 2 ( A« , 2 1 H SIDE CUN 

1 DUCT OR RINGS,5X)/8X,5HTEMP.,7X,6(1H(, A2, 1H) , 1 1 X ) , /7X , 6hDEG. K, 

2 2C15H CONDUCTIVITY) ,2(30H CONDUCTIVITY RESISTIVITY)/ 

3 1 i X , 2 ( 1 5H (WaTT/CM/K) ) * 2 (28H (WATT/CM/K) (OhM*CM) , 2X ) ) 

WRITE (6,253) ( T TT ( I ) , ZK IN ( I ) , ?KSS ( I ) , ZK CRH ( I ) , RCRH ( I ) , ZKCRC(I) » 

1 RCRC (I), 1*1, 15) 

RETURN 

END 
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SUBROUTINE COUPLE 

COMMON /MDCPQP/*A0(9) ,0R(8) , TEMP ( 9 , 300 ) , VOC (300) , RPC ( 300) , PE (300) , 

1 QT (300) , NOPTIM, VPC(300) ,RLPC(300) »DTM0D(300) , CUR (300) , LNRAD (9) * 

2 ZLP, ZLN,ZLI,ZlNLP,ZLPNI,ZLTE,NGTi,NGT2» I TPERT , PCMULT , NODUMP, IVD, 

3 QGEN,aGENL,QTE,Q0UT,ZlD,Z0D,ZKEND,ZK9,ZKR,RPPC,RPN,T0LTEM,TR£j, 

4 VREQ,P£REQ, ZLREQ, WTCON, ZLPN, IOUMZ , RN, OGRF , DGTRF , DGIRF , R AOK , JDUMP, 

5 DP(30),NC 

COMMON /M0DCPL/CURM0D(15) , I T ( 300) , DT TE , CONRN, CONRP, ZL IP, ZNC , 

1 ALRb4,MAXTEM,QGAM,CGAM, I NDI , R INT 4 , R INTb , RC0N4 , RCONb, C 1 FEH, C 1FEC , 

2 CONST 1»C0NST2,CQNST3, CONS T4, CONST 5, NR50PT»SEGTE, I CRMDL»RAOC 
COMMON /DATIN/TTT(15),ZKMIC1 (15) ,ZKMIC2(15) , ZKSS ( 1 5) , ZKBN ( 1 b) , 

1 ZKIN(15),Z2NA(15),Z2NR(15),22NK(15),Z2PA(15),Z2PR(15),Z2PK(15J, 

2Z2NRK(l5) ,Z2PRK(15),Z3NA(15) ,Z3NR(15),Z3NK (15) ,Z3NRK(15) ,Z3PA(l5) , 

3 Z3PR(i5) ,Z3PK(15) ,Z3PRK (15) #ZKCRH(15) ,2KCRC(15) ,RCRH(15) *RCRC(15) 
4,RHP(15) »RPRH0 

COMMON /TITLE/FETM(3) ,CASCA0(8,2) #SSIN(3) ,CMT (25) , COND (3, 5) ', TEMKF » 
1 C«M0(2,2) , I3N, I2N, I3P, I2P 
COMMON /CPLRIT/DT(8#300) # OPTZLP(300),OPTCPL(300) 

DIMENSION XZK(9) f ZK(9) ,TC0N(11),ZKEQ(8) ,B(100) 

REAL LNRAD 


OATA XZK/9*0,/,LNDUMP/0/ f 

ZKTEP/O./, 

zkten/o,/ 

, T50P/670 

,/,ZKEQ/8*0./ 

DATA B/ 6HTFEH 

, 6HZKP 1 

,6HRINT4 

,6HZKDT 

, 6HTFEC 

,6HZKP2 , 

16HRINT6 ,6HQSAVE 

,6HF£RH 

, 6HZKN1 

, 6WTCL0H 

, 6HQGEN 

, 6HF ERC , 

26HZKN2 , 6HTCLDC 

, 6HQGAM 

,6HALP1 

, 6HZK 1 1 

, 6HTBNH 

, 6HQTQMP 1 » 

36HALP2 , 6HZK I 2 

, 6HTBNC 

, 6HQT0MP2 

,6HALN1 

,6HZKT£P 

, feHDTMOD , 

46HQT0MN1 ,6HALN2 

,6HZKT£N 

, 6HDTTE 

,6HQT0MN2 

, 6HZKR 

, fcHZKTE , 

56HT0LTEM,6H0T0M 

, 6HRH0P 1 

,6HZKTEQ 

, 6HTC0N 

, 6HQC 

,6HRH0P2 , 

66H2KECPC,6HALT1 

, 6HQP 

, 6HRH0N1 

,6HZLP 

, 6HALT2 

, 6HQ J , 

76HRH0N2 , 6HZLN 

, 6HALNH 1 

, 6HQ IN 

, 6HRHQP 

,6HZLT£ 

, 6HALNH2 , 

86HQ0UT , 6HRH0N 

, 6HRP 

,6HALNC1 

,6HVPC 

, 6HC0NST 1 

,6HRN , 

96HALNC2 ,6HRPC 

,6HC0NST2,6HRFEH 

, 6HALPH1 

, 6HV0C 

, 6HC0NST3, 

A6HRFEC , 6HALPH2 

, 6HCUR 

, 6HC0NST 4 

, 6HR I NTFH 

, 6HALPC 1 

,6HPE , 

16HZKH0T , 6HRI NTF C 

, 6HALPC2 

» 6HQTE 

,6HZ*CDLD 

, 6HRPPC 

,6HRLPC , 

26HQT , 6HQ6QL 

, 6HQB0L 

, 6HTREDL 

,6HTRB0L 

, 6HPCMULT / 


T0LTEM.AMAX1 (TOLTEM, ,001) 

nseg«segte 

YNCiZNC 

IF (NGT1.GE.3) YNC*PCMULT 
ZKECPC«ZKENO/YNC 
C0NST6.0,5*ZLR£Q/(ZLI*YNC)-l .0 
DO 6010 1*1, NC 
QSAVE*0. 

ITN0W*0 
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C INITIALIZE radial temperatures 
IF (1N0I.6T.1) GO TO 4053 
IF (I.fO.l) GO TO 4045 
11*1-1 

IF (QGEN.LE.O,) GO TO 5005 
DTMODCl)*DTMODCU) 

TEMPU, I)«TEMP(9,I)-*-0TM00(I) 

5005 DTMnDR»OTMOD(l)/OTMOO(Il) 

DT(1»I)*DT(1»I1) aDTMQDR 

do aoai j s i , 7 

J9*J+1 

TEMPtJq, IJaTEMPCj, D-DT (J, I) 

4041 DT(J9,I)5DT(J9,U)*DTM0DR 
GO TO 4040 

4045 IF (OGEN.LE.O.) GO TU 5004 

DTMUD(l)«0.5*GGEN*C0NRN/SI<TTT,Z2NKr TEMP(9, I) » 15#IND# 1, 1) 
TEMPtl, I)*TFMP(9,I)+DTM00(I) 

5004 TEMP (5, n*TEMP(l,I)-LNRAQ(4)*DTM0D(I)/AlR64 
DU 4051 J*2,4 
TEMPI J, n«TEMP(l,I) 

4051 TEMPIJ+4, I)*TEMP(9, I) 

DO 4052 J* 1 » 8 

4052 DT t J, I)*TEMP(J f I) -TE MP C J+ 1 , I ) 

GO TO 4040 

4053 IF (NGTl.LT. J.UR.ITPERT.GT.l) GO TO 4040 
IF (ITPERT.EQ.OJ GO TO 4055 

lTPBKTmO 

DTM0D(l)*TEMP(l#n-TEMP(9,n 

DTM0DReDTM0D(l)/DTM0D(3) 

DU 8013 J t ■ 1 » NC 

DO 4054 J* 1 » 7 

DT (J, J1)*DT(J, Jl)*DTMODR 

4054 TEMPI J+l , J1)*TEMP( J, J1 )-DTCJr Jl) 

6013 0T(8, Jl )*0TM0PR*0T(8, Jl) 

4040 DTTE*TFMPC4,I)-TEMPCb,I) 

IF (PTTE.GT ,0TM0D (I) ) DTTE * DTMOD(I) 

IF (NR50PT.EQ.0) GO TO 4060 

TEMP (5, I)«AMIN1 (AMAX1 (T5UP, TEMP (4, I) ), TEMP (5,1)) 

DT (4, I)«TEMP(4, I) -TEMP (5, I) 

DT (5, I)*TEMP(5» I ) -TEMP (6, I) 

4060 0ELT1«DT (4, I) /SEGTE 
DELT2*DT(5,I)/SEGTE 
TT1 s TEMP (5,1) ♦ , 5*DEL T 1 
TT2 s TEMP (6, I) ♦ ,5*0ELT 2 
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Astronuclear 

Laboratory 


ALP 1 * 0. 

ALP? = 0, 

RKP1 • 0. 

RKNl * 0. 

RKP2 a 0. 

RKN2 b 0, 

ZKP1 s 0, 

ZKP2 a 0. 

ALN1 s 0. 

ALN2 = 0, 

ZKNl a o, 

ZKN2 * 0, 
ZKIUO. 

ZK I2«0, 

DO 20 jal » NSEG 


ALPl 

s 

ALPl 

♦ 

SICTTT, 

Z3PA, 

TT1, 

15, IND, 1,1) 

ALP 2 

3 

ALP2 

♦ 

SICTTT, 

Z2PA, 

TT2, 

15,IN0,1,1) 

RKP1 

O 

RKP 1 

♦ 

SICTTT, 

Z3PRK 

, TT 1 

,15, INO, 1,1) 

RKP2 

B 

RKP2 

♦ 

SICTTT, 

Z2PRK 

,TT2 

, 15, IND, 1 » 1 ) 

ZKP1 

S 

ZKP1 


SICTTT, 

Z3PK 

,TTt 

,15,lND,l,l) 

ZKP2 

8 

ZKP2 

+ 

SICTTT, 

Z2PK 

, TT2 

, 15, IND, 1, 1) 

ALN1 

8 

ALN1 

+ 

SI CTTT , 

Z3NA 

a TT 1 

,15, IND, 1,1) 

ALN2 

B 

ALN2 

♦ 

SICTTT, 

Z2N A 

, TT2 

,15, IND, 1,1) 

RKN1 

S 

RK N 1 

♦ 

SICTTT, 

Z3NRK 

a TT 1 

,15,IND,1,1) 

RKN2 

s 

RKN2 

♦ 

SI CTTT, 

Z2NRK 

, TT2 

,15, IND, 1,1) 

ZKNl 

8 

ZKNl 

♦ 

SICTTT, 

Z3NK 

, TT l 

,15, IND, l, 1) 

ZKN2 

8 

ZKN2 

♦ 

SICTTT, 

Z2NK 

,TT2 

,15, IND, 1,1) 


ZKIlaZKlltSI(TTT,ZKMICl,m,15,IND,l,l) 
ZKI2sZKl2+SICTTT,ZKMIC2»TT2,15,IN0,l,l) 
TT1 a jTl + OELTl 
20 TT2 a yT2 ♦ 0ELT2 

alpi*alpi/segte 

ALP2»ALP2/SEGTE 

alnualni/segte 

ALN2aALN2/SEGTE 

ZKPUZKP1/SEGTE 

ZKP2aZKP2/SEGTE 

ZKNIpZkNI/SEGTE 

ZKN2bZKN2/SEGTE 

RKPIbRrPI/SEGTE 

RKP2«RKP2/SEGTE 

rkni«rkni/segte 

RKN2aRKN2/SEGTE 

ZKllaZKll/SEGTE 
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2kI2*ZkI2/SEGTE 

RHUPlsRKPl/ZKPl 


RH0P2*RKP2/ZKP2 

RHONlsRKNi/ZKNl 

RH0N2srKN2/ZKN2 

«055 ALPH1 s SHTTT,Z3PA,TEMP(4,I),15,IND,1,1) 

ALNHl s SICTTT»Z3NA»TEMP(4,I),15fIND'l'i) 

ALPH 2 s SI(TTT,Z2PA,TEMPf5,n»l5,IN0,l,l) 

ALNH2 8 SI(TTT,Z2NA,TEMPC5,I),15,IN0,1,1) 

ALPC1 = $I (TTT,Z3PA, TEMP ( 5, I) , 15, IND, 1, 1) 

ALNCI*SICTTT, Z 3* A, TEMP (5, I), 15, IND, 1,1) 

ALPC2*Sl (TTT,Z2PA,TEMPC6, I) , 15, IND, 1 ,1) 
ALNC2 ssKTTT,Z2NA,TEMP(6, I ) , 1 5, IND, 1 , 1 ) 

UPsTEHPU, 1)* CaLPH1 + ALNH1)+TEMP(5» I)*(ALPH2+ALNH2*ALPC1«ALNC1) 


CALCULATE OPTIMUM SEGMENTING RADIUS, RAO (5) , IF SPECIFIED 
IF CNR50PT.EU.0J GO TO 4061 
RAD (5) bRAD (6) 

IF CDT(5,IJ.LE.O.) GO TO 4062 
TT1*ZKn 1*DT (4,I)/fZKN2*DT(5,I)) 

TT2=1 .O/CTTl+l.O) 

WAD(5)8(RAD(6)**CTT1*T T2))*CRAD(A)**TT2) 

IF (HAD(5) .LT.RA0C4) .0R.RADC5) .GT.RADCbJ ) GO TO 4063 
4062 LNRAQ(4)*AL0G(RAD(5) /RAD (4) J 
LNKA0(5)*AL0G(RAD(6)/RA0(5) ) 

4061 WHOP * CINRAD(5J*RH0P2 ♦ INRAD (4) *RHOP 1 J /ALR64 
RHON * CLNRAD (5J *RH0N2 + LNRAD (4) *RH0N1 ) / ALR64 
ZKTEP s ALR64 /(LNRaD(5)/ZKP2 ♦ LNR AD ( 4 )/ IK P 1) 

ZKTEN a ALR64/(LNRAD(5)/ZKN2 ♦ LNRAD ( 4} /ZKNl ) 
CALCULATE OPTIMUM ZLP IF NO VALUE HAS BEEN READ INTO Z(12) 
GO TU (4042,4001,4002), IDUMZ 
4001 TT2sKHQN*ZKTEP/(RH0P*ZKTEN} 


IF (TT2.GE.0,) GO TO 1096 
4063 JNDIbO 

N0DUMPs5 
GO TU 3112 

1096 ZLNLP5S0RTCTT2) 

4002 IF (IVd.NE.2) GO TO 1097 

ZKI«ALR64/(LNRA0(5) /ZK 12 + LNRaD (4) /ZK I 1 ) 

ZLNLP s ZLNLP*SQRT((cONST6-ZKI/ZkTEN) /(C 0NST6-ZKI/ZKTEP) ) 
ZLN«ZLl*C0NST6/(,5+,5/ZLNLP) 

1097 ZLPsZLn/ZLNLP 

zlip=2.o*zli/zlp 

CQNST5sl5,9593*ZLP 
ZLPN« ZLN+ZLP 
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Astronuclear 

Laboratory 


ZLPNIsZLPN^ZLl 

ZLTEfZLPNI+ZlI 

QGEN*0g£NL*ZLT£ 

C0NST2s15.9593*ZLTE 

OPTZLP(I)»ZLP 

GPTCPL(I)*ZLTE 

0PTCPL(I+100J8YNC*ZLTE 

RINT4*RC0N4*ZIPN/ZLP 

RINT6*RC0N6*ZLPN/ZLP 

CONRPicONSTl/ZUP 

IF (ICRMDL ,NE. 0) GO TO 4005 

ClFEH*l5.9593*0R(3)*(RAD(3)+RA0(4))/ZLTE 

C1FEC*15.9593*0R(6)*(RA0(6)*RAD(7))/ZLTE 

GO TO 4042 

CAUCULATE RADIUS OF CIRCULAR CONDUCTOR RING ELECTRICAL STREAMLINE 
4005 RCH0T« A MIN1 (DR(3) ,ZLN,ZLP) 

RCOLD-AMINl (DR(6),ZLN,ZLP) 

C1FEH»5.08*RAD(4)*ALOG(3, 14159*RCH0T/ZLU1 .0) 
C1FEC=5,08*RAD(6)*AL0G(3.14159*RC0LD/ZLI+1.0) 

CALCULATE COUPLE RESISTANCE 
4042 RINTFH 8 RINT4*TEMp(4, I)*TEMP(4, I) 

RINTFC«RINT6*TEMP(6# I)*TEMP(6, I) 

TFEC « TEMP (7,1) ♦ ,5*0T(6,I) 

TFEH « TEMP (4,1) ♦ ,5*DT(3,I) 

FERC * SI(TTT,RCRC,TFEC,15,IND,1,1) 

FERH a SI(TTT,RCRH,TFEH, 15, I NO, 1 , 1 ) 

RFEHsFERH/CIFEH 
RFEC8FERC/C1FEC 
RN « C0NRN8RH0N 
RP s CONRP*RHOP 
RPN8RP+RN 

TTE«,5*(TEMP(a, I)+TEMP(6, I) ) 

RPPC8RpRH 0*SI(TTT,RHP,TTE» 15, IND, 1 , 1 ) / YNc 

UJiRFEH+RlNTFH+,5*RPN 

RPC(I) 8 QJ+,5*RPN+RFEC+RINTFC+RPPC 

CALCULATION DP THERMAL CONDUCTANCES 
TBNH8 TEMP (3# 13 ♦.5*DT(2,I) 

TBNC* TEMP (8# I ) ♦,5*DT(7,I) 

TCLDH * TEMP (2, I ) +,5*0T(1,I) 

TCLDC » TEMP(9,I) +.5*07(8,1) 

XZK(l) s SI(TTT,ZKIN,TCLDH,15,IND,1,1) 
XZK(8)*SI(TTT,ZKSS,TCL0C, 15, I NO, 1, 1) 

XZK ( 2)8 SI (TTT,ZKBN,TBNH, 15, I NO, 1 , 1) 

XZK (7) b SI(TTT,ZKBN,TBNC, 15, INO, 1, 1) 
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XZKC3)a(SI(TTT,ZKCRH,TFEH, 15 ,InD,1,1)*ZLPNI+SI(TTT,ZKMiC1,TFEH, 

1 15,InD,1,1)*ZLI)/ZLTE 

XZK(6) 5 (SI(TTT,ZKCKC,TFEC, 15, InO,1,1)*ZLPNI*SI(TTT,ZKMIC2,TFEC, 

1 15, InD, 1, l)*ZLI)/ZLTE 

ZK(l)acONST2*XZK(l)/LNRADCl) 

ZK(2)acUNST2#XZK(2)/(LNRAO(2)*(l,0+XZK(2)*C0NST3)) 

2K(3)>CON3T2*XZK(3)/LNRA0(3') 

ZK (4) 8 C 0NST5*(ZLNLP*ZKN1 + ZKP1+ZI.IP*ZKH ) /L NR AD (4) 

IF (LNR AD ( 4 ) .LE.O.) ZK(4)*1,0F30 

ZK (5) SC0NST5* (ZLNLP*ZKN2+ZKP2*ZHP*ZKI2I/LNRAD(5) 

IF (LNRAD(S) ,LE.O.) ZK(5)«1.0E30 
ZKTE * ZK(4)*ZK(5)/(ZK(4)+ZK(5)) 

ZK(6)a C 0NST2*XZK(6) /LNRA0C6) 

ZK(7)ac0NST2*XZK(7)/(LNRAD(7)*(l,0+XZK(7)*C0NST4) ) 

ZK (8) s C 0NST2*XZK (8) /LNRA0C6) 

IF (ZKr.GT.O.) Zk(8)sZK(8)*Zk9*ZlTE/(ZK(8)*ZK9*ZlTE) 

AUU * ALPl + AtNl 
ALT2 a ALP2+ALN2 

CALCULATION of voltages and current 

VUC(I) ■ DTf«,n*ALTi+DTC5,I)*ALT2 

IT(I)aiT(I)M 

ITNUW=tTNOrt+l 

ITCONaiTNOw-MAXrFH 

ITCONaMAXOdTCfJN, 1) 

ocsdtte*zkte 

QENDaZKECPC*DTMOD(n 

GO TO (4047,4048,4049,4044,4050), NRT2 

4047 CURCDsO. 

GO TO 4046 

4049 RLPC ( I) aHPC (I) 

GO TO 4048 

4050 RLPC (I) aHPC (1) aSQRT (1 , + VOC (I) * (OP-»VOC f I ) *QJ/RPC (I) ) / (RPC ( I ) * (QC+ 
1 QEND))) 

GO TO 4048 

4044 QTEsPCMULT*VOC(I)-VREQ 
RLPCdm ,OE30 

IF (QTE,GT.O.) RLPCmsVRFOaRPc(I)/OTE 

4048 CUR(l)aCURMOD(INOI) 

IF (NGTl ,GE,3,0W,INDI,E0. 1 ) CUR (I ) a VUC ( I ) / (RPC ( I ) +RLPC ( I ) ) 

4046 VPC(I)aVOC(I)-CUR(I)*RPC(I) 

CALCULATE ENERGY TERMS 

QT0MP1.CUR(I)*(ALPH1*TEMP(4, I)-ALPC1*TEMP(5, I)-ALP1*DT(4, I ) ) 
UT0MN1«CUR(I)*(ALNH1*7EMP(4,J).ALNC1*TEMP(5,I)-ALN1*DT(4,I)) 
QT0MP2«CUR(I)*(ALPH2*TEMP(5, I)-ALPC2*TEMP(6, I)*ALP2*DT (5, I) ) 
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QT0MN2aCUR(I)*(AtNH2*TtMP(5, I)-ALNC2*TEMP(b» I)-ALN2*0T(5, I) ) 

QT0M3QT0HN1+0T0MN2+QTQMP1+QT0HP2 

IF CQTOM.GT.O.) QT0M*,5*GT0M 

QPsCURf 1)*QP 

QJ*CUR(I)*CUR(I)*QJ 

QTEsQP*GJ-QTOM*QGAM*CGAM 

QIN»GC«fOTE 

PE(I) a VPC(I)*CUR(I) 

CALCULATE EFFICIENCY 
QT ( I ) «Ql N+GGAM 
GOUT«OT(I)-PE(I) 

2KH0T«l,0/Cl,/ZKcn + l./ZK(2)+l./ZK(3)) 

ZKC0LDsl./Cl./ZK(6)4l,/ZK C7)*1./ZK(B)) 

IF (GGEN.GT,0,.AN0.VREQ,GT.0..aN0.0GRF.EQ,0,) GINb,5*(QGEN*0END 
1 4-gin) 

CALCULATE TEMPERATURE DROPS AND NEW RADIAL TEHp, PROFILE 

zkdt*qin/qc 

2KTEQ«*ZKTE*ZKDT 

ZKEQ(1)«ZK(1) 

ZKEG (2 ) »ZK (2) 

ZKEQ ( 3) »ZK (3) 

ZKEG(4)«ZK(4)*ZKDT 

ZKEQ(5)»ZK(5)*ZKDT 

ZKOIsQlN/GOUT 

ZKEQ(6)3ZK(6)*ZKQT 

ZKEQ(7)«ZK(7)*ZKDT 

ZKEG(8)«ZK(8)*ZKDT 

IF (ZKR.GT.O.) TeMP(9,I)b(ZKR*TREJ^00UT+ZKECPC*TEHPC1,I))/ 

1 (ZKR+ZKECPC) 

IF CRA0C.LE.0..0R,INDI.EG.1.0R.N0PTIM,LT t 4) GO TO 5008 
5008 IF (GGEN.GT.O.) GO TO 5001 

ZKDTwDTMODCI) / (1 ,0/ZKHOT +1 ,0/ZKTEQ+QOUT/ (QIN*ZKCOLD) ) 
DT(l#I)aZKDT/ZK£0(l) 

DO 500fe J«2,8 
J 1 • J * 1 

TEHP(J f I)»TEHP(Ji # I)-DT(J1»I) 

5006 DT(J#I)*ZKDT/ZKEQ(J) 

TCON ( ITCON) «ABS ( 1 ,-GSAVE/QT CD) 

GSAVESQT(I) 

GO TO 5007 
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CALCULATE MOT CLAD TEMPERATURE IF HE AT GENERATION RATE IS SPECIFIED 
5001 ZKDTsQgEN-QEND 
00 500? J = l,8 
J9*9»J 

DT(J9, l)aZK0T/2K£0(J9) 

500? TEMPCJq, I)«TF.MP(J9+t,l)*6T(J9, I) 

TCUN(ITC0N)«? # *AhS(1 ,-QSaVE/TEhP(1* I) ) 

QSAVEstEMPU, IJ 
OTMODCl)aOSAVF-TEMP(9, I) 

5007 IF CN0DUMP.E0.0.0R.N0DUMP.EQ.2) GO TO 3007 
3112 WRITE (6/1 00) (CMT(J), J«l,25), INDI,IT(I), I 
100 FORMAT (1H1,4X,2A4, IX, 18 a4,2H ( , F2 , 0 , 2 ( 1 H/ , F2 , 0 ) , 2H) ,2A4// 

1 20X/19H CURRENT HERAT ION, 15, 23h, TEMPERATURE ITERATION, 

1 15, UH, ON COUPLE, I5//1 IX, 3HRAD, 1 IX, 4HTEMP, 12X,2HDT, 12X,3HXZK, 

2 13X,2hZk, 1?X,4HZKEQ, 13X,5HLNRaD) 

WRITE (6,3002) (J.RAO(J) , TEMP ( j , I ) , DT ( J , I ) , XZK ( J ) , ZK ( J) , ZKEQ ( J ) , 

1 LNRAO(J) , j*i.a) 

3002 FORMAT CI3,1H., 7G15.5) 

J s 9 

WRITE (6,3003) J , RAD ( J ) , TEMP ( J , I) , OTMOD (I) , ZK9, PCMULT, ALR64 

3003 FORMAT (13, 1H.,3G15.5, 15X , 3G15.5//) 

WRITE (6,3004) 0(1) , TFEH , p (2) , 7KP 1 , B ( 3 ) , R I NT4 , B (4) , ZKDT 
WRITE (6,3004) H(5),TFEC,B(6),ZKP2,B(7),RINT6,6(8),0SAVE 
WRHE (6,3004) B(9) ,FERH,B(10) ,ZKNl,B(U) ,TCL0H,8(12) ,l)GEN 
WRITE (6,3004) 8(13) ,FERC,B(14) ,ZKN2,B(15) ,TCLDC,B(16) , OGAM 
WRITE (6,3004) 8(17),ALPl,B(18),ZKIl,B(l9),T8NH,B(20),QT0MPl 
WRITE (6,3004) B(21),ALP2,B(22),ZKI2,B(2i),TBNC,B(24),oTOMP2 
WRITE (6,3004) B (25) , ALN 1 , B (26) , ZK TEP , B (27 ) , OTMOD (I ) , B ( 28) , QT OMN l 
WRITE (6,3004) 8 ( 29 ) , ALN2 , 8 ( 30 ) , ZKTEN , H ( 3 1 ) , DTTE , B ( 32 ) , 0T0MN2 
WRITE (6,3004) B(33),ZKR ,B(34) , ZKTE ,8 (35) , TOl TEM, B (36) ,UTOM 

WRITE (6,3004) B(37) ,RH0P1,B(38) , ZKTEQ, B ( 39) , TCON ( ITCON) , B (40) , QC 
WRITE (6,3004) 8(41),RH0P2,B(42),ZKECPC»B(43) ,ALT1,B(44) ,0P 
WRITE (6,3004) B(45),RH0Nl,B(46)f ZLP,B(47),ALT2,B(48),0J 
WRITE (6,3004) 9 (49 ) , RH0N2 , B (50 ) , ZLN, B (M ) , ALNH1 , B (52) , Q IN 
WRITE (6,3004) 9 (53 ) , RHOP , 8 ( 54 ) , ZL TE , B ($5) , ALNH2, B ( 56 ) , QOUT 
WRITE (6,3004) B (57 ) , RHQN , B (58) , RP , B (59) , ALNC 1 , B (60) , V PC ( I ) 

WRITE (6,3004) 8 (6 l ) , CONST l , 8 (62) » «N , B (63 ) , ALNC2 , B (64 ) , RPC ( I ) 
WRITE (6,3004) B(65),C0NST2,9(66),RFEH,B(67), ALPH1 ,B(68) , VOC(I) 
WRITE (6,3004) 8(69) ,CONST 3,8(70) ,RFEC,B (71) ,ALPH2,B (7?) ,CUR( I) 
WRITE (6,3004) 8 ( 73 ) , C0NST4 , 8 (7 4 ) , R I NTFH , B (75) , ALPC 1 , B ( 76 ) , PE ( I ) 
WRITE (6,3004) B(77) ,ZKH0T,B(78) ,RINTFC»B(79) , Al PC2,B(B0) ,QTE 
WRITE (6,3004) B ( 8 1 ) , ZKCC1L D , B ( 82) , RPPC , B ( 83 ) , RL.PC ( 1 ) , B ( 84 ) , UT ( I ) 

3004 FORMAT (2X,4(A6* 1Hs,G12,5, 10X) ) 

IF (INDI.EO.O) WRITE (6,3006) 
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3006 FORMAT ( 1 HO , 2 1 X * 66HTH 1 5 CASE DUMPED* THEN SUPPRESSED BECAUSE OPTlM 
1UM ZLP IS IMAGINARY) 

IF (NODUMP.EQ.l. AND, ITCON.LT.il) GO TO 3007 
WRITE (6*3008) (TCON(J) *J»1, ITcON) 

3008 FORMAT (10H0 TCON * * 10 (F9,6* 1H, ) ) 

GO TO 9000 

3007 IF CTCQN(ITCON) .GT.TOLTEM) GO TO 3009 
IF (NODUMP, EQ. 2) GO TO 3112 

GO TO 4000 

3009 IF (ITCON, LE, 10) GO TO 4040 
GO TO 3112 

4000 IF (N0PTIM.NE.5,0R,I,NE, 1) GO TO 6010 
C * SET UP B.O.L, PARAMETERS FOR CURIUM MISSION CALCULATIONS 
IDUMZ* 1 

QGEN*1,211*(QT(1)'MJEND) 

IF (INOI.EQ.l) GO TO 6010 
N G T 2 * 4 

QEOLsPCMULT*(QOUT«-QEND) 

QBQLsPcMULT* (OGEN-PE (2) ) 

ZKDf ■G bOL/UEOL 
TRE0L*TEMP(9, D-RADK 

TRBOL»(ZKDT*TREOL**4-(ZKOT*1,)*TREJ)**,25 

TEMP(9,3)«TEMP(9,2) 

TF.MP(9,2)»TRB0L + ZKDT*RA0K 
ZKDT*TEMP(9,2)*TEMP(9,3) 

DO 6012 J«l*8 

6012 TEMP(J,2)pTEMP(J,2)4-ZK0T 

IF (NODUMP. EQ.O) GO TO 6011 

WRITE (6,3004) 8 (85 ) , QEOL * B ( 86 ) * 080L * B (87 ) , TREOL , 0 C 88 ) , TRBOL 
WRITE (6*3004) B(89)*PCMULT 
6011 CALL PHONY(DP* J0UMP,0,) 

6010 CONTINUE 
RETURN 
END 
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BLOCK DATA 

COMMON /DATIN/TTf (15),2KMIC1 (15),ZKMIC2(15), ZKSS U 5 ) , ZKBN ( 1 b ) , 

1 ZK IN ( 1 5) » Z2NA(15) , Z2NR(15), Z?NK (15) ,Z2PA (15) ,72PK(15) , Z2PK (1 5) , 
2Z2NRK(i5) ,Z2PRK(15),Z3NA(15),Z3NRC15) , Z3nK ( l 5) , Z 3NRK ( 1 5) , Z3PAC15) , 
3 Z3PRU5) »Z3PK(l5)f23PRK(15')»ZKCRH(15)>ZKCRC(15)rRCRH()5) ,RCRC(15) 
4, «HP(15) #RPRH(1 

COMMON /F0DAT/ALPHA(15,B),RH0(15,8),ZKUN(16,8),ZKSSIN( 16, 3) , 
l RFE TM ( 1 5, 4), ZKFETM Ub, 3 ),ZKMICA( 16) 

COMMON /TITLE/FETM(3),CA8CAD(8,2),SSIN(3),CMT(25),C0N0(3,5) ,TPMKF, 
1 CRMD(?,2), I3N,l2N,nP,I2P 

C TEMPERATURE AT WHICH EACH PARAMETER IS EVALUATED (DEG.K) 

DATA tTT /300., 350. ,400. ,450. ,500. ,550. ,600. ,650. ,700., 

1 7 50., 800., 850,, 900. ,950,, 1000,/ 

C SEEBECK COEFFICIENT OF TEGS-3N MATERIAL ( VOL TS/DEG . K ) 

DATA aiPHA/ 8.526E-05, 1 .0S8E-04, 1 .254E-04, 1 .41 IE-04, l .656E-04, 

2 1 .803E-04, 1 .960E-04,2,097E-04,2.205E*Oa,2,254E-04, 

3 2.254E*04,2.185E-04,2,009E-04, 1 .735E-04, 1.303E-04, 

C SEEBECK COEFFICIENT OF TEGS-2N AND GE-NL MATERIAL ( VOL TS/DEG.K ) 

1 1.150E*04, 1.369E-04, 1 ,5B7E-0a , 1 , 806E-04 , 2 . 005E-04 , 

1 2.1 T6E-04,2.338E*04,2.452E-0«,2.509E-04,2.«90E-04, 

1 2.386E-04,2,205E-04, 1.977E-04, 1.730F*04, 1.464E-04, 

C SEEBECK COEFFICIENT OF TEGS-3P MATERIAL ( VOLT S/DEG ,K ) 

1 5.0l5E*05,5.807E-05, 8.974E-05, l . 1 40F -04 , 1 , 3 94E -Oa , 

1 l.b68E-04, 1.91 lE-04,2.122E-04,2.281E-04,2.«07E-04, 

1 2.502F-04,2.555E-04,2.544E-04,2,4bOF-04,2.?81E-04, 

C SEEBECK COEFFICIENT OF TEGS-2P MaTfRIAL ( V0LTS/PEG , K ) 

1 4,940E-05,7,980E-05, 1.102E-04, 1 . 387 E -04 , 1 . 748E-04 , 

1 2.090E-04,2,375E*04,2.594E-04,2,698E-04,2,726E-04, 

1 2.70eE-04,2,660E-04,2.594E*0«,2,508E*04,2.394E*04, 

C SEEBECK COEFFICIENT OF RCA-NB MATERIAL (VOlTS/DEG , K ) 

1 1 .030E-04, 1 .210E-04, 1 , 4 1 OE-04 # 1 , 6 1 0E-04 , 1 ,B2oE-04, 

2 2.020E-04, 2.230E-04, 2.450E-04, 2.b00E-04,2.bb0E-04, 

. 3 2,b20E-04,2.530E-04,2,430E-04,2,300E-04,2.14oF-04, 

C SEEBECK COEFFICIENT OF RC A GE-TE P-TYPE MATERIAL (VOLTS/DEG. K) 

1 0,52E-04,0.78E-04, 1 .05E-04, 1.35E-04, 1 , 7 3E-0« , 2 . 05E-04 , 

1 2.38E-04,2 t 6ftE-04,2,83E-04,2 t 98E-0«,3.1 IE-04, 

1 3,09E-04,3.0QE-04,2.19E*04, .48E-4/ 

C RESISTIVITY OF TEGS-3N MATERIAL (OhM-CM) 

DATA RhO / .0002628, ,0003623, ,0004967, ,0006722, ,0008992, 

1 .0011785, .0015277, .0019206, .0023658, .0028198, 

1 .0032476, .0035880, .0037714, ,0037452, .0036142, 

C RESISTIVITY OF TEGS-2N AND GE-NL MATERIAL (OhM-CM) 

1 .0004132, ,0005932, ,0008393, .0011702, .0015737, 
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1 


,0020822, 

,0026552, 

.0032685, 

,0038738, 

,0043984, 

1 


.0046728, 

,0047615, 

.0046082, 

.0042209, 

,0037124, 

RESISTIVITY 

OF 

TEGS-3P MATERIAL (OhM 

-CM) 



1 


.0009940, 

.0010799, 

.0012395, 

.0014972, 

.0018408, 

1 


.0022090, 

.0026753, 

.0031539, 

.0036939, 

,0043811 , 

1 


,0051542, 

.0057310, 

.0061360, 

.0063814, 

,0064428, 

RESISTIVITY 

OF 

TEGS-2P MATERIAL (OhM 

-CM) 



1 


.0004031, 

.0005493, 

.0007445, 

,0009987, 

,0013165, 

1 


.0017342, 

.0022517, 

.0028237, 

.0035227, 

.0042400, 

1 


.0049482, 

.0055656, 

.0059923, 

,0061739, 

,0059469, 


c resistivity of RCAcNB material (ohm-cm.) 

1 4.750E-04,6.300E-04,8.250E-04, 1 .070E-03, 1 .390E-03, 

2 1.820E-03,2,360E-03,3.040E-03,3.800E-03,4.570E-03, 

3 5.300£-03,6.000E-03,6.450E-03,6,800E-03,6,750E-03, 

C RESISTIVITY OF RCA G£-TE P-TYPE MATERIAL (0HM*CM) 

1 3.7E-04, 4,8E-04,6,SE-04,9.7E-04, 13.8E-04, 19.4E-04, 

1 25.5E-04,32.8E-04,40.5E-04, 48 , 5E-04 , 58. 0E-04 , 67 , OE-04 , 

1 77.56-04, 8b, 0E-04, 98. 0E-04/ 

C THERMAL CONDUCTIVITY OF TEGS-3N MATERIAL (W AT TS/CM ,/DEG .K ) 

DATA Zk ON / ,0928220, ,037551b# .0332035, .0289872, .0250344, 

1 ,0213451, .0181829, .0lb20b5, .0148889, .015020b, 

1 ,0tb20S5, ,0178558, .0191052, .0210816, .023058, .298, 

C THERMAL CONDUCTIVITY OF TEGS- 2 N AND GE-NL MATERIAL t W AT TS/cM/DEG . K ) 

1 ,0286942, .0244433, .0208299, .0180888, .0157287, 

1 ,0139220, .0127530, .012221b, .0128593, ,014134b, 

1 .0159412, ,0178542, ,0204098, . 0233805 ,. 026037 4 ,, 298 , 

C THERMAL CONDUCTIVITY OF TEGS-3P MATERIAL [WATTS/CM, /DEG. K) 

1 ,0375000, .0305000, .0248000, .0199000, .0180000, 

1 .0130000, .0114000, .0110000, .0112000, .0119000, 

1 .0135000, .0186000, .0208000, .0265000, . 035000 ,. 298 , 

C THERMAL CONDUCTIVITY OF TEGS-2P MATERIAL ( Wa TTS/CM. /DEG . K ) 

1 .0428878, ,0348821, .0281394, ,0227591, .0182988, 

1 .0148678, .0130379, .0125804, ,0128092, .0136097, 

1 .0154396, .0189850, .0235597, .0303074, .0400286, .298, 

c thermal conductivity of rca-n(a) material (watts/cm, /deg.kj 

1 .0200, .0179, .0161,, 0143,. 01 27,. 01 11,. 0101,. 00971, .00974, ,0101 , 

2 ,0106, ,0113, ,0122, .0133, .0147, ,296, 

C THERMAL CONDUCTIVITY OF RCA-P(A) MATERIAL (WATTS/CM, /DEG.K ) 

1 .0345, ,0270, .0222, .0189, ,01 56, ,01 35, .01 18, .01 08, ,01 03, .0103, 

2 .0109, .0128, .0181, .0303, .0653, .298/ 

C THERMAL CONDUCTIVITY OF 80R0N NITRIDE ( WATT /CM/DEG .K ) 

DATA ZKBN /. 303, .301, .299, .296, .294, .291, .289, .286, .284 

l , .282, .279, .277, .274, .272, .270, ,081/ 

C THERMAL CONDUCTIVITY OF IRON ( w A T T/CM/DEG. K ) 
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data ZKFETM /. 725, .696, .664, .635, .608, ,57a, .540, .510, .479 

1 , .449, .420, .390, .360, .335, .313, .280, 

C THERMAL CONDUCTIVITY OF TUNGSTEN f w A T T/CM/QE g . K ) - TPP C DATA 

a 1,78,1.70,1.62,1.56,1.51,1 .46, 1,4 1,1.3 t, 1.34 

3 , 1 .31, 1 .29, 1.26, 1.24, 1.22, 1 .20, ,700, 

C THERMAL conductivity of MOLYBDENUM (WATT/CM/DEG.K) • TPRC data 

4 1.3 7,1.35, 1.33,1.31,1.30,1.27,1.26, 1.24,1.22 

5 , 1.20, 1.19, 1.17, 1.15, 1.14, 1.13, .368/ 
c resistivity UF IRON (OHM*CM) 

DATA RFETM /l.l5E*5,1.40E-5,l.73E-5,2.08E-5,?.52E-5, 

1 2.97F-5,3.52F>5,4.12E-5,4.886-»5,5.60E-5,6.34E-5,7.l0E-5, 

2 7.91E-5,8.80E-*5,9.82E-5, 

C RESISTIVITY OF TUNGSTEN (OHM*CM) TPRC DATA 

3 5.00E-6,6.25E-6,7.50E-6,8.70E*6, 1 .00E-5, 1 , 1 2E-5, 1 .25E*5, 1.38 

4 E-5, 1 .53E-5, 1 .67E-5, 1 .82E-5 , 1 , 9BF-S, 2 , 1 5E*5 , 2 . 29E-5, 2 . 43E-5 , 

C RESISTIVITY OF MOLYBDENUM (OHM*CM) TPRC DATA 

5 5.50E-6,6.50E"6,7,SOE»6,8.50E-6,9.50E-6, 1 . 04E-5 , 1 . 1 3E*5 , 1 , 25 

6 E-5, 1 .37E-5, 1 .47E-5, 1 .58E-5, 1 .69E-5, 1 .BOE-5, 1 .90E-5, 2.00E-5, 

C RESISTIVITY OF NICKEL (OHM*CM) • 

7 .000011, .000013, .000015, ,0000173, .0000 1 98 , , 0000226 , . 0000257 , 

8 .000029, .0000319, .0000346, . 000037 , .0000392, , 00004 1 , , 0000427 , 

9 .000044/ 

c thermal conductivity of stainless steel (watt/cm/deg.k j 

DATA /KSSIN / . 1 45 , , 1 57 , . 1 85, . 1 7 2 , , 1 7 9 , . 1 84 , . 1 92 , . 1 98 , . 205 

1 , .212, .219, .225, .232, .239, .246, .28 3, 

C THERMAL CONDUCTIVITY of INCONNEL (wATT/CM/OEG.K) 

3 .117, .126, .135, .143, .152, .161, .170, .179, .188 

1 , . 197, .205, .215, .223, .230, ,240, .296, 

C THERMAL CONDUCTIVITY of T-lll ( T P R c PUWDUF UNIV.) 

1 .418, .431 , .439, .452, .464, .481 , .490, .502, .515, .527, .540, .552, 

2 ,565, .577, .586, .600/ 

C THERMAL conductivity of MICA (NATT/CM/DEG.KJ 

data ZkMICA /, 0290, .0310, , 0330, . 0342, .0336, .0306, .0282, .0268, .0262 
1 , .0262, .0268, . 0268 , . 026 8 , . 0 268 , .0268, .105/ 

DATA CASCAD/4HT£GS,4H GE »2*4HTEGS,2*4HRCA-,2*1H ,4H 3n ,4H NL , 

1 4H 3P ,4H 2P , 4h N , 4H P ,2*1H / , C0ND/4H0PE N , 4H CIR,4HCUIT, 

2 4H FIX,4HED L,4H0AD , 4HM ATC , 4HHED , 4HL 0 AD , 4HF I XE , 4HD V-,4HL0AD, 

3 4H0PTI , 4HMUM , 4 HL HAD/ , F E TM / 2HFE , 1 HW , 2HMq / , SS I N/2HSS , 2h I N , 2HT A / , 

4 CMT/25*4H**W,CRMD/4H C IR , 4WnUlT, 4H M0 d,4HULE / 

END 
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SUBROUTINE OPTIM ( ISECT , I RTRN) 

DIMENSION VPL(4#5) ,0MHD(4) ,DRCP(2»2) ,PDWT(2,4) ,VPL1 (3) ,PWCT(4) 

1 ,VRL(5)#DRCP1 (2),SKP(12) 

COMMON /MODOPT/IPDWT, DR4 , DR INcR » I R I TE , DRTE , JOPTIM, I Z9 , I ENQK , IPUn , 
1 RIlNC,R0lNC,RlMAX,R9MAX,DTMIN f QREQ, THING, TCI NC,TCM AX, THMAX 
2,NPERT,ZNP£RT,CKPP,CKFT#CTE0 C#CTEIC#DELTl»RIPERT,RLMAX,ZNMIN,Z9SaV 
COMMON /TITLE/F6TM(3) #CASCAD(8,2) »SSIN(3) ,CMT (25 ) , COND ( 3 , 5) , TEMKp , 
I CRMD (2#2) , I3N, I2N , I3P, I?P 

COMMON /M0CP0P/RAD(9) ,D«(a) , TE mP( 9, 300} # VOC (300) , RPC (300) r PE (300) , 

1 QT(300) »N0PTIM» VPC(300) ,RLPC(300) »DTM0D(300) ,CUR(300) ,LNRAD(9)» 

2 ZLP,Z|_N,ZLI,ZLNLP,ZLPNI,ZLTE,NGT1,NGT2* itpert,pcmult, NODUMP, I VD, 

3 OGEN,QGENL#OTEf Q0UT,ZID,Z0 D#ZkEND#ZK9,ZkR,RPPC#RPN,T0lTEM,TREJ# 

4 VREQ,PERE0,ZLREQ, wtcon, zlpn, IDUMZ,RN,DGRF,DGTRF,DGLRF,RADK, jdump, 

5 DP(30),NC 

DATA VPL/9H V , 4H0C ,4H ( VO, 4HLTS) , 4H CUR» 4HRENT , 4H (AM,4HPS) , 

1 4HCKT.,4HLEN.,4H (I»4HN,) ,4h RL»4H0AD »4H (0H,4HMS) , 4H RL» 

2 4H0AD ,4H(M-0,4HHMS)/,PoWT/4HPDEN,4HW/CC#4H WT , , 4HLBS , , 4H CUR,4hA 
3MPS , 4h T5 ,4HDEG,/,0MHD/3H M«,1H , 1 HC# 1 HH/ , DRCP/4H DR5, 4HDRTE , 

4 «H N0.,4HCPIS/,SKP/1H ,1H0,8*1H 
EQUIVALENCE (VPLl(l),OCV), (VPL1(3)#ZCIR) 

REAL LNRAD 

IF (ISECT, EQ,2) GO TO 3051 

C - INITIALIZE VARIABLES USED IN OPTlMZATION ROUTINE 
IRTRNsi 

IF (JOPTIM.NE.O) GO TO 4014 
JOPTlMsl 
IPERTso 
PCPERTsNPERT 
PCSAV*PCMULT 
RLOADaPCMULT*RLPC(i) 

IRITEa^ 

IV0«1» TEMPERATURE PARAMETRIC WITH SPECIFIED GEOMETRY 
IVD«2, TEMPERATURE PARAMETRIC WITH SPECIFIED LENGTH ANr) VOLTAGE 
IVD«3, TEMP, PARAMETRIC WITH SPECIFIED LENGTH AND NO, OF COUPLES 
I VD*4, PARAMETRIC ON NO, COUPLES WITH SPECIFIED VLOAD AND LENGTH 
I VD®5» PARAMETRIC ON NO, COUPLES WITH SPECIFIED VLOAD AND POWER 
klom*iooo, 

IDRCP«2 
IPDWT»1 
ISKPi2 

ZCIR«PCMULT*ZLTE 
C - INITIALIZE PARAMETERS FOR CURIUM STUDY CALCULATIONS 
N0P28N0PTIM/2 
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IF (N0P2.NE.2) GO TO 6000 
OC v so « 

Q s 0 . 

IF (NUPTIM.EQ.b) NC*2 

1 V PL a 3 

IVRLsl 

TREJaTREJ*TR£J*TR£J*TREJ 
OGTRf S V ,/OGTRF 
TEMP(1,2)«TEMP(1.1) 

TEMP(9,2I*T£MP(9* I ; 

ISKP*2 
1 PO W T s ^ 

GO TO 4009 

C - INITIALIZE PARAMETERS for LOAD CURVE CALCULATIONS 
6000 IF (RLMAX.LE.O.) GO TO 4000 
IRMLsl 
ISKPsl 

RLSAVSRLPCU) 

I V PL s 1 
1 VWLsS 

IF (RLOAD.LT.i*) GO TO 3999 

RLOMsl .0 

1VRL=4 

3999 IPL)W T s 3 
IRLP=l 

IF (WLPERT.LE.O.) RLPERT aRLSA V 

IF ( (RLMAX-RLSAVJ/RLPERT.GT.200.) GO TO 4009 

IRLPsO 

GO Tu 4009 

4000 IF (ZLREQ.LE.0. I GO TO 3998 
ZLNa.5*ZLREQ/PCMULT-ZLl 
ZLP»ZLN 

1 VD»3 
ZCIRsZlTE 

3998 IF (VRfU.LE.O,) GO TO 4002 
IVDs2 

RLPCUJsO. 

IF (NPfRT.LE.O) GO TO 4002 
NGT1 *4 
NGT2s4 
I V 0 a 4 
I PDWT s 3 

RLPC(l)*l.OE10 

IF (PEREG. LE.O.. AND. 0RE9. LE.O.) GO TO 4002 
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I VD*5 
PCPERTs 1 , 

CUREQsPEREQ/VREQ 
4002 I VRU * 1 

IF (NGT2.GE.4) I VRL*5 
I VPL * 1 

IF (IVRL.EQ.l) I VPL*2 
I R5PR T *>2 

IF (DRINCR.LE.O.) GO TO 4004 

IF (N0PTIM.NE.2) IR5PRT«1 

IPDWTstt 

IDRCPsl 

GO TO 4009 

4004 IF (N0PTIM.NE.3) GO TO 4009 

C - INITIALIZE VARIABLES FOR TEMPERATURE DERIVATIVE CALCULATIONS 
TCINC«13,89 
THINCsTCINC 
thmaxbtempu, n ♦thing 
TCMAX aTEMP(9 f D4-TCINC 
TEMP (1, !)«TEMP(1,1)-THINC 
TEMP (9, 1)* TEMP (9,1) -TC INC 

4009 T1H0T«TEMP(1, 1) 

T ICOLDaTEMP (9,1) 

IF (IPUN.EO.l) WRITE (7,101) (CMT ( J) , J*3, 23) 

101 FORMAT (17A4,A2# 1X,F2.0,2(1H/,F2.0)) 

J*38 

IF (THMAX.EQ.O) GO TO 4038 
I8(THMAX«T1H0T)/THINC+1 , 

IF (N0PTIM.GT.3) I »2* I 
J*(39/I)*I 

4038 LCNTM«AMAX0(2, J) 

ZLNSAVsZLN 

ZLPSAVsZLP 

C - INCREMENT 0.0, IF INITIAL VALUE VIOLATES MINIMUM T/E THICKNESS 
4014 IF (ISECT.EQ.3) GO TO 4023 

WRITE (6,100) (CMT(J),J«1,25),CASCAD(I3N,2),CASCAD(I2N,2), 

1 CASCAD(I3P,2),CASCAD(I2P,2) 

100 FORMAT (1H1,4X,2A4,1X,18A4,2H ( , F2.0, 2 ( 1 H/ , F2 . 0) , 2H) ,2A4/ 

1 47X,1H(,A3,1H*»A3,1H,,A3,1H-,A3,1H)) 

WRITE (6, 4010) ZI0,Z00»DR(1) ,0RTE»0R(8) 

4010 FORMAT (13X,4HI.0.,8X, 1Hb,F 8.4,4H IN, , 45X , 4H0 ,D . , 8x » 1 H» , F8 , 4 , 

14H IN./26H INNER CLAO THICKNESS » , F8 * 4 , 4H IN. , 9X , 6HQRTE «,F7.4, 

2 4H IN., 10X,22H0UTER CLAD THICKNESS e F8.4,4H IN. ) 

IF (ZCIR.LE.O.) GO TO 4061 
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ZLPNI*ZtP*ZLN + ZU 
WTaWATE (P) 

J=1 

IF (wTcON.GT.O. ) J = 2 

IF (NGT2.NE,2,0R.RLPERT,GT.0.) GO TO 4040 
1=2 

IF (RUlAD.GE.l.) GO TO 4045 

rloao*ioqo.*rlqad 
1 = 1 

4045 WRITE (6,4041) ZC I R, RLOAD , OHHD ( I ) , CRMD ( 1 , J ) , C RmD ( 2 , J ) , U T 

4041 FORMAT C 1 OX , V 6HC I«CU I T LENGTH «,FB,3,4H IN. , 6X , 7HRL0AU *,F9.3, 

1 A3, 4H0HMS, 1 1X,2A4,9H HEIGHT «,F8.3,5H IBS, ) 

GO TO 4042 

4040 WRITE (6,4043) ZC IR , (CONO ( I , NGT2) , I « 1 , 3) , CRMO ( t , J ) , CRMO (2 » J ) , W T 

4043 FORMAT (10X, 16HCIRCUIT LENGTH m,F8.4,4H IN . , 1 1 X , 3 A4 , 1 8X , 2 A4 , 

1 9H WFXGHT =,F8*3»5H LBS.) 

GO TO 4042 

4061 IVPL=3 

WRITE (6,4063) (COND ( I , NGT2 ) , I a 1 , 3) 

4063 FORMAT (49X.3A4) 

4042 WRITE (6,4044) ORCP ( 1 > I DRCP ) # f VPL ( I , I VPL ) , I * 1 , 2) , ( V PL (I , I VRU # 

J 1=1,2) , (CRM0C1, IENDK) , I ■ 1,2), PONT (1, IPOWT) , TEMKF , TFMKF , D«CP (2, 

2 I DRCP ) , C VPL Cl* I VPL), 1*3,4), (VPL(I, JVRL) , I =3 , 4) , PDWT ( ? , I PQW T) 

4044 FORMAT ( 1 7 HO rHOT TCOLO, 4X , A«, 1 4H ZLN ZLP , 4X , 2 A 4 , 2X , 

1 2A4,4X,30HRGEN VLOAD P-OUT 0, A 4 , 5X , 3HE T A , A 2 , 6X , A 4 / 

2 2 ( 7H (DEG.,A1,1H)),3X,A4,18H (IN.) (IN.) ,2A4,2X,2A4, 

3 54H (M-OHMS) (VOLTS) (WATTS) (WATTS) (PCT.) (,A4,1HJ) 

LCNTsO 

I ZLP= 1 

ITHlNCsl 

ITCINC«0 

etamaxbO, 

return 

3051 IF (N0P2.NE.2) GO TO 6001 
I TPER Ts2 

C - ITERATE ON NO, OF COUPLES AND ZLN FOR CURIUM MISSION CALCULATIONS 
I0UMZ=2 
NG T2*5 
QGENsO. 

IRTRNa? 

P=PCMULT*PE(1) 

VLD=PCMULT*VPC(1 ) 

IF ((ABS(l,-0CV/VLD)+ABS(l.-Q/p)).LT..()l,0R.IPERT.GT.10) GOTO 4001 
IPERT=IPERT+1 
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PCMULTsIFIX(VREQ/VPC(1))+1 

ZLNe(PEREQ*ZlPN/(P£(n*PCMUl,T*Cl.M./ZLNi.P)) + ZLN)*.5 

ocv*vld 

QaP 

IF (N0PTIM.EQ.4J TEMP(l,2)aT6Mpn#l)+20. 
0a0GLRF*EXP(166l6.*(0GTRF*l ./TEHP(1#2)) )/ZLN 
OaAMINl (D, ,75) 

0P(12)sD*(RPC(NC)+0P(12)*(RPC(NC)"RN))/(RN*(l.-D)) 

IF (N0PTIM.EQ.5) DP(12)sO, 

CALL PHONY (OP»JDUMP f 100000.) 

RETURN 

6001 QCV*PCmULT*V0C(1) 

IF (IVD.EQ.l) GO TO 4001 
IF (IZLP.NE.2) GO TO 4035 
toltem«io.*toltem 
GO TO «001 

4035 TOLTEM.O. 1*T0LTEM 
IZLP*2 
IRTRN*2 

IF (IVr>.EQ.2) GO TO 4028 

4027 IF (PEREQ.LE.O.) GO TO 4018 

IF (OCV.LE.VREQ) GO TO 4029 

ZLNaZLN/(l.+(VOC(l)-VREQ/PCMULT-CUREQ*RPC(l) ) / (CUREQ*RPN) ) 
if (Zln.gt.o.) Return 
PCMULTsPCMULT+1 
GO TO 4037 

4018 IF (UReQ.GT.O.) ZLRE9sCQREU-ZKEND*DTM 0D(1))*ZLTE/QT(1) 
ZLNa(ZLREQ/PCMULT-2.0*ZLI)/(l , 0 +l.O/ZLNLP) 

4020 IF (VREQ.LE.O.) RETURN ' 

CALCULATE MINIMUM NO, OF COUPLES TO ACHIEVE SPECIFIED LOAD VOLTAGE 
PsOCV^vREQ 

IF (P.LE.O.) GO TO 4029 
RLPC(1)8VRE9*RPC(1)*ZLP*ZLNLP/(P*ZLN) 

return 

40129 PCMULTbIFIX(VREQ/V0C(1)+2,0) 

toltem«io,*toltem 
I ZLP* 1 

RPPC»RPPC/PCMULT 
RPC ( 1 ) aRPC ( 1 ) -RPPC*PCMULT+RPPC 
4037 QCVaPCMULT*VOC(l) 

PE(l)avREQ* (OCV^VREQ)/ (PCMULT*PCMULT*RPC(1) ) 
RLPC(l)«VREQ/(PCMULT*CUREQ) 

GO TO 4027 

CALCULATE NO. OF COUPLES TO ACHIEVE. SPECIFIED LOAD VOLTAGE 
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4026 PCMULT=IFIX(VRE0/VPC(1)+1.0) 

ZLN« (ZL«tQ/PCMUtT-2.0*ZLI)7 (1 .0 + 1 ,0/ZLNLP) 

PET URN 

4001 TEMPHFxDK2FK(TEMp(l,lJ) 

TEMPCFs0K2FK (TEMP (9, 1 } ) 

IF CIZlP.EQ.2.0«.0ELTL.EQ.0.) GO TO 4064 
CALCULATE total LENGTH TO ACHIEVE specified clad MISMATCH 
IZLPsP 
IRTRNsp 

ZC1K = DELTL/CCTEIC*(TEMPHF^70,)-CTE0C*(TEMPCF-70.) ) 

IF (ZLREQ.GT.O. J ZCIR»ZLREU 

PCMULT*IFIX(ZCIR/(ZLNSAV*(1,0M,0/ZLNLP)+2,0*ZLI))M 

ZLNs(ZcIR/PCHULT-2.0*ZLn / (1 ,0+1 ,0/ZLNLP) 

return 

4064 ETAUaPE{l)/rQT(l)+ZKENO*OTMOO(i)/PCMUlTl 
IF (FTaO.LT.ETAMaX) GO TO 4015 

zcih*zl,te*pcmult 
RLUADapCMULT*RLPC (1 ) 

VRL(IVrL)=PLOH*RLOAD 

VRLtnsPCHULT*VQc(l) 

DRCP1 (1 1 a 100. *DR (5) /DRTE 
DRCP1 (2)*PCMULT 

V0L1 * 51.48*ZCIR*RAD(1)*RAD(1) 

Vl,DaPCMULT*VPC(l) 

UHMsPCmULToRPC(I) 

PsPCHUlT*PEC1) 

U*PCMUi.T*QT(1) 

VPL1 C2)»CUR(l) 

IF (1P0WT.E9.2) PWCT CIPOwT) sWATECQMOD) 

PWCT (3) «CUR{ 1 ) 

PWCT [4) S0K2FK (TEMP (5, 1 ) ) 

IF (N0PTIM.NE.2J GO TO 4003 
ETAMAXaETAO 

QMODsZlP 

PWCT(1)*ZLN 
GO TO 4017 
4015 ZLPaQMfiO 

ZLNaPWnT(l) 

ETAOaETAMAX 

4003 QMOOsQ + ZKEND*OTMOO (1 ) 

PWCT(l)sUMOO/VOLl 

IF ( IPUN.E0.3) WRITE (/,40bb) TEMPHF, TE«PCF, QMDD, OCV, VlD, CUR ( 1 ) # 0, 
1 (CMt(j), J»21,2i) 

4055 FORMAT (7HREF MOD , 1 2X , 2F7 . lr F8 . 2 , 3F8 . 4, F8 . 2, 3 A2/ 1 H ) 
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IF (IPUN.EG.l) WRITE (7 , 4050} Z ID# 200# TEMPHF, TEMPCF , ZC IR , VLD, OHM, 
1 ET AO# P 

4050 FORMAT (lX,2F7.3,2F7 t l,5E10.3) 

IF (N0PTIM.NE.3) GO TO 4031 
JOPTIMsJOPTIM+1 

V0C(J0PTIM)*0CV/(DTM00(1)*CKFT) 

RPC ( JOPTIM) *OHM 
QT (J0PTIM)bDTM0D(1)*CKFT/QM0D 
4031 lcnt*lcnt*iskp 

IF (LCnT.LE.LCNTM) GO TO 4033 

WRITE (6,100) (CMT(J),J*1,25),CASCAD(I3N,2),CASCAD(I2N,2), 

1 CASCAD(I3P,2) ,CaSCAD(I2P,2) 

WR I f E (6 , 40 1 0) ZlD,ZQD,DR(l),DRTE,DR(a) 

IF (Rl-OAD.IT . 1 « ) GO TO 4034 
RL0M*1 .0 
I VKL * 4 

VRL ( 4 ) aRLO AD 

4034 WRITE (6,4044) ORCP ( 1 , IDRCP ) , ( VPL ( I , I VPL ) , I a 1 , 2 ) , ( VPL ( I r I VRL ) , 

1 1=1,2) , (CRMD ( 1 , IENDK ) , I * 1 , 2 ) , PDWT ( 1 , IPDWT ) , TEMKF , TEMKF » ORCP (2, 

2 IDRCP), (VPL (I, I VPL), I =3, 4) , (VPLCIrlVRL) , 1*3,4), PDWT (2,1 PDWT) 
PCPERTs2.0*PCPERT 

4033 WRITE (6,4011) SKP( I SKP ), TEMPHF , TEMPCF , DRCP 1 ( IDRCP) , ZLN , ZLP , 

1 VPL1 (I VPL) , VRL (I VRL) ,OHM, VLO , P , QMOD , ET AO , PWCT ( IPDWT ) 

4011 FORMAT (Al,F7.1,F9.1,F7.0,2F8.4,2Fi0.3,3pF10,3,2(0PF10.3), 

1 FlO, 1,2PF10,4,OPF10.3) 

LCNT«2 

GO TO (6005, 4016, 4025, 6006, 6002), NOPTIM 
C - WRITE B.O.L. PARAMETERS FOR CURIUM MISSION CALCULATIONS 

6002 WRITE (7,4050) Z ID , ZOO , TEMPHF , TEMPCF , ZC I R , VLD , PWC T (2 ) , F T AO , P 
TEMPHFsDK2FK(TEMP(l ,2) ) 

TEMPCF«DK2FK(TEMP(9,2)) 

0CV*PCMULT*V0C(2) 

0HM=1000.*PCMULT*RPC(2) 

vld*pcmult*vpc(2) 

P*PCMUlT*PE(2) 

QM0D*PCMULT*QT(2)+ZKEND*DTM0D(2) 

ETAO«P/QMOD 

Q*QM00/V0L1 

WRITE (6,6003) SKP(IPERT) , TEMPhF , TEMPCF , 0, 0, PWCT (2) , OCV, OHM , VLD, 

1 P,OMOo,ETAO, IPERT 

6003 FORMAT (IX, Al,Fll.l,F9.1,2PF6.1,0PF8.2,F9.2,Fl0.3,F9.1,F10.2, 

6 F10.1,F11.1,2PF9.3,I6) 

LCNT=LcNT*l 

CALL Ph0NY(DP,J0UMP,100000.) 
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IPEK1 =0 
GO TO 4055 

6006 WRITE (6,6006) 0, PWCT ( 1 ) , PWCT (2) , IPERT 
6006 FORMAT (22X,2PF6.1,0PF8.2,F9.2,59X,Ib) 

IPERT *0 
GO TO 4025 

6005 IF (RLMAX.LE.O.) GO TO 4017 
C - INCREMENT LOAD RESISTANCE, IF SPECIFIED ( Z ( 47 ) , GT , 0 , ) 
RLPC(i)aRLPC(l)+RLPERT 
IRTRN*5 
I T H I NC * 1 

GO TU (5001,5002,5003), I RML 

5001 IF (RLPC(t),UT.RPC(l)) GO TO 5003 
IRMEa2 

RSAV = RlPC (1 ) 

N G T 2 * 3 
RETURN 

5002 I RMC * 3 
N G T 2 » 2 

RUPCClJaRSAV 

5003 P=R.99*RIPERT 

IF (IRlP.EO.UANO.RLPC(I) .GE.P) RLPERT*10.*RLPERT 

IF (RLPC(l).LE.RLMAX) return 

WLPC(l)aRLSAV 

LCNTsSo 

RLUMal 000. 

IVRL*5 

IRMLsI 

if (IRLP.EQ.l) RLPERTsRLSAV 
GO TU 4025 

4017 IF (DRINCR.LE.O.) GO TO 4021 
C - INCREMENT OR (4) IF SPECIFIED (NZ (32) .GT.O) 
[)R(4)eoR(4)+ ORlNCR*ORTE 
RA0(5)*RA0(4)+DR(4) 

DR(5)»RAD(6)-RA0(5) 

IRTRNS3 

IF (OR (4) .GT.DRTE) GO TO 4008 

LNRAD(4)=AL0GCRAD(5)/RAD(4)) 

LNRA0(5)*AL0G(RA0(6)/RA0(5)) 

RE TURN 

4008 IF (N0PTIM.EQ.2) GO TO 4003 
WRITE (6,4011) SKP (2 ) 

LCNTsLCNT+1 
4016 OR (4) sqR4 
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RAD(5)sRAD(4)4-DR(4) 

DRC5)*RAD(6)-RA0(5) 

ETAMAXsO. 

4021 IZ L P * 1 

IF CNPERT.LE.O.J go TO 4025 

C - INCREMENT AXIAL DIMENSIONS, IF SPECIFIED CNZ(2b) .GT.O) 

IF (PCSAV.GT.l) GO TO 4051 

C - INCREMENT NUMBER OF COUPLES IF PcMULT (NZC12)) .LE.l) 

IF (IPUN.EQ.2) WRITE (7,817) (C«T (J ), J*3 , 1 5) , PCMULT , PCMULT , ZLN , ZLP 
817 FORMAT (13A4, 1H,F4.0,8H C0UPLES/6H 1 ll2,F6,0/6H 2 111,2G12.5) 
PCMULTsPCMULT+PCPERT 

IF (ZLN.GE.ZNMIN.AND.PEREQ.LE.O.) GO TO 4026 
IPERT * IPERT+ 1 

IF (CUREQ.GT *0, ) RLPC (1) ■VREQ/ (PCMULT*CUREQ) 

IF (IVD.GT.4.AN0.IPERT.LE.NPERT) GO TO 4036 
IPERTsO 

PCMUUTaPCSAV 

WRITE (6,4011) SkP(2) 

GO TO 4025 

4026 Psl,/(l.*QTE*ZLTE*RLPC(l)*PCPERT/(RPC(l)*PCMULT*(ZLTE*t)T(l)-2,*ZlI 
1*(QT(1 J-QTE))) J 
ZLN*P*ZLN 
4036 IZLP-l 
IRTRN*2 
return 

C - INCREMENT ZLN, IF SPECIFIED 

4051 IF (IPERT. LE.NPERT) GO TO 4052 
IPERT *0 

ZLN*ZLNSAV 
ZLP*ZLPSAV 
WRITE (6,4011) 

GO TO 4025 

4052 IPERT»IPERT+1 
ZLN«ZLN+ZNPERT 
ZLP»ZLP-ZNPERT 

IF (ZLP.LE.O.) GO TO 4023 

IRTRN«2 

RETURN 

4025 ITPERTnl 

DTM0D(3)aDTM0D(l) 

IRTRN*3 

IF (THINC.LE.O.) GO TO 6004 
C - INCREMENT HOT CLAO TEMPERATURE, IF SPECIFIED 
6009 TEMP(1, 1)«TEMP(1, D+THINC 
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ISKP=1 

IF (TEMP(1, 1) .LE.THMAX) RETURN 

TEMP(l, DaTIHOT 

ISKP=2 

ZCIRsO. 

faOOa IF (TCINC.LE.O.) GO TU 4023 
C - INCREMENT COLO clad temperature, if specified 

TEMP (9, 1)*TEMP(R, D+TCINC 
IF (TEmP(9, 1) .GT.TCMAX) GO TO 6007 
IF ((TEMP(1,1)-TEMP(9,1)).GT.0TMIN) RETURN 
IF ( (ThMaX-TEHP( 9# l) ) .GT.DTMIN) GO TU 6009 
600/ TEMP(9 , I)sTJCOLD 
ZCIRsO. 

C - INCREMENT outer t/f washer THICKNESS 
9023 IF (TEmPCI, 1) .gT.THMAX) GO TO 3062 
IF (IZ9.E0.1) GO TO 4019 
OR(5)=nR(5)+R0INc 
IF (DR(5) .GT.R9MAX) GO TO 3062 
IRTRN*4 
RETURN. 

C - increment outer diameter 

4019 RAD(9)«RAD(9)+R0INC 

IF ( R At> ( 9 ) .GT.RRMAX) GO TO 3062 

I R I R N * 6 

RETURN 

C - INCREMENT INNER DIAMETER 
3062 RADCl)*RAD(l)+RlINC 

IF (RADC1) .GT.RlMAX) GO TU 4024 

I RTRNs 7 

RETURN 

4024 IF (NOP T IM, NE, 3) GO TO 4030 
CALCULATE ANo WRITE TEMPERATURE DERIVATIVES, IF SPECIFIED 
Ps27.78*CKFT 

V0C(2)«(V0C(7)-V0C(5))/P 
V0C(4)s(V0C(9)-V0C(3))/P 
RPC (2)* (RPC (7) -RPC (5) ) /P 
RPC (4) s(RPC(9)-RPC (3) )/P 
aT(2)«(QT(7)*QT(5))/P 
QT(4 )s(QT(9)-QT(3))/P 
WRITE (6,4074) 

4074 FORMAT (1H0/36X,34HCALCULATED TEMPERATURE DERIVATIVES) 

WRITE (6,4071) 0MHD(4) , V0C(2) , TEMKF, TEMKF 
4071 F0RMAT(1H033X,6HPAL/*T,A1,2H »,E13.«,12H VOLTS/DEG, , A 1 , 5H/0EG. A 1 ) 
WRITE (6,4071) 0MH0(3) ,V0C(4) , TEMKF, TFMKF 


D-78 




Astronuclear 

laboratory 


WRITE (6, 4072) OMHD (4) , RPC (2) , TEMKF 

4072 F0RMAT(1H033X,6H*RG/#T,A1,2H «,E13.4,11H OHMS/DEG . t A 1 ) 

WRITE (6, 4072) OMHD (3), RPC (4) , TeMKF 

WRITE (6,4073) OMHD ( 4) , QT (2) , TEMKF , TEMKF 

4073 FORMAT (IH033X,6H#T I /»T, A1,2H «,E13,4,6H DEG. , Al , iOH/WATT/OEG . , A 1) 
WRITE (6,4073) OMHD (3) , QT ( 4) , TEMKF , TEMKF 

TEMPO, 1)«TEMP(1,1)+THINC 
TEMPO, 1)*TEMP(9, D+TCINC 
4030 IRTRNsl 
RETURN 
END 
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SUBROUTINE LIFE(OP) 

COMMON /MDCPOP/RAOCR) ,DR(«) ,TEMP(9,300) r V0C(300) ,RPC(300) ,PE(300) , 

1 QT(30o) ,NQPTIM, VPC(300) , RLPC ( 300) , DTMOO C 300) , CUR C300 ) , l NR AO ( 9) , 

2 ZLP,ZLN, ZLI, ZLNLP#ZLPNl,ZLTE,NGTl,NGTPr ITPERT»PCMUI T , nODUMP , I VD , 

3 QGEN,QGENL,QTErQOUT,ZID,ZOD, Z kEND,ZK 9, ZKR, RPPC , RPN, TOl TEM, TRE J • 

a VREU,PEREQ, ZLREO, WTCON,ZI_PN, If)UMZ,RN, DGRF, PGTRF, DGLRF, RADK , JDUMP, 
5 0P(3O) f NC 

COMMON /MODCPL/CURMOO(15) , I T ( 300 ) , DTTE , CUNRN , CONRP , ZL I P » ZNC , 

1 ALK64, MAXTEM,QGAM,CGQM* INOI,RlNT4,RlNT6,RCON4,RCONb,ClFEH,CiFEC, 

2 CONST 1 , CONST? ,C0NST3,C0NST4, CONST*, NR50PT, SEGTE, IC'RMDL#R a DC 
COMMON /MODOPT/lPDWT, DR4.DRINCR# IR I T£ , OR TE , JOPT IM , I Z9 , I ENOK , IPUN, 

1 RIINC,R0INC,R1MAX, R9MAX,D T MIN, OREO, THING, TCI NC»TCM AX, THM AX 
2,NPERT, ZNPERT,CKFP,CKFT,CTEnc,cTElC,DELTL#RLPERT,RLMAX,ZNMIN,Z9SAV 
COMMON /TITLE/FETMC3) ,CASCA0(8,2) ,SSIN(3) ,CMT(25) ,C0ND(3,5) , TEMKF, 
1 CRMD(2,2) , I3N, I2N, I3P, I2P 

COMMON /molif/hRs,hrinc,hhmax,hflf»radpct, ITCR, ITHQ 
DIMENSION CIRM00I2) 

DATA C iRMOD/aHC IR , A H M [1 D / 

QEPC=ZkEnD * DTMOD(l ) /pcmult 

if c joptim.ne.oj go to no 

JQP T I M a i 
THSAVstEMP( 1 , n 
TCSAV = TEMP(9, 1) 

WRITE (6,100) (CMT(J) , J=1 ,25) 

100 FORMAT C1H1 , 4X,2A4, IX, 18A«,2H ( , F2 . 0 , 2 ( 1 H/ , F2 , 0 ) , ?H ) ,?A4) 

WRITE (6,101) CIRMOD(IENDK),CIRMOD(IENDK) , TEMKF, TEMKF 

101 FORMAT (//92rt0 TIME THOT TCOLO R1.0AD RGEN 

1 VOC VLOAD CURRENT POWER 0-»A«,7H ETA-,A1/ 

2 9X,bHHOURS,2(0H DEG., A 1 ) ,2(9H M-0HMS),2(9H VOLTS), 

3 9H AMPS , 2 ( 9H WATTS), 9h PCT./) 

UBOLsO. 

IF (IThO,EQ,1.AND,HFUF.EU.O.) GO TD 102 

QBUlsOT(l)+GEPC 

QGENsOBOL 

102 IF (ITCR.NE.I) ZkRb(Q0UT+0EPC)/(TEMP(9, 1)-TREJ) 

110 TEMPHFsDK2FK (TEMP(1 , 1 ) ) 

TEMPCFsDK2FK(TEMP(9,1)) 

RLUADbPCMULT*RLPC(1) 

UHM=PCMULT*RPC(1) 

C *C UK ( 1 ) 

OCVsPCMULT*VOC(l) 

VUD = PCMULT*VPC (1 ) 

PePCMUL.T*PE (1 ) 
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Q?PCMULT*(QT(1)+QEPC) 

e*p/q 

WRITE (6, 111) HRS.TEHPHF, TEHPCFf RLOAD»OHM,QCV, VLD,CUR(1) , P,Q,E 
1 ITU) 

111 FORMAT (1X,F13,0,2F9.1,3P2F9.3,0P4F9,3,F9.1,2PF9.3,I10) 
JOPTIMbJOPTIM+1 

PE(J0PTIM)«P 

V0CCJ0PTIM)*HR3 

hrs*hrs+hrinc 

IF (HRS.GT.HRMAX) go TO 112 
IF (HFlF.GT.O.) QGENbOBOU*£XPC«HRS/HFLF) 

CALL PHONY(DPr JDUMP,HRS) 

RETURN 

112 J0PTlM*-2 

IF (RAqPCT.EQ.O,) RETURN 
Q*l,0Efe*(l.-PE(3)/PE(2))/V0C(3) 

Esi.OEb*(i.«PE(7)/PE{2))/VOC(7) 

WRITE (6,201) V0C(3),Q,V0C(7),E 
201 FORMAT ( 1 H0//25X , F 1 0, 0, 25H HOUR POWER DEGRADATION «,F8,3, 
l 22H PC T # PER 10,000 HOURS) 

TEMP ( 1 , 1 ) *T HSA V 
TEMP (9, 1)«TCSAV 
RETURN 
END 


\ 
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SUBROUTINE PUMPCIRTRN) 

COMMON /MOOpPT/lPOWT, ORA, OR INC R, J«ITE,0RTE» JOPTIM, IZ9, IENOK, IPUN, 
1 RIINC,ROINC, RiMAX,R9MAX,0TMIN,QRE9»THlNC,TCINC»TCMAX r THMAX 

NPERT, ZNPeRT,CKFP,CKFT,CTE0C»CTElC*0ELn.»RLPE«T,RLMAX,ZNMIN,Z9SAV 
COMMUN /MDCP0P/RA0(9) ,DR(8 ) , TEmP(9,300) , VOC ( 300 ) , RPC ( 3o0 ) » PE ( 300 ) , 

1 01 (300) ,N0PTIM,VPC(300) , RLPC ( 300) , DTMUP ( 300) , CUR (300) , LNRAD ( 9) # 

2 ZLP,ZLN,ZLI,ZLNLP,ZLPNI,ZLTE,NGT1,NGT2# I TPERT , PCMUL T , NODUMP , IVD, 

3 0GEN,UGENL#QTE*Q0UT # Z1D. ZOD, ZkEND,ZK 9, ZKR,RPPCf RPN,'rOLTEM*TRFJ» 

a VRF 0 ,PERE 0 #ZLRe 0 »WTC 0 N,ZI,PM, lDUMZ,RN,OGRF,OGTRF,DGLRFrRAOK, JDUMP, 
5 DP ( 30 ) , NC 

COMMON /TITLE/FETM(3) , C A SC AD (8, 2) , SSIN(3) ,CMT(25) ,C0N0(3,5) , TEMKF, 
1 CRMD(2,2) , I3N, I2N, I3P,I2P 
C0MMUN/TEPMP/NPUMP,DR8»C0D2 
IF (JOPTIM-1) 102, 103,201 
C - INITIALIZE VARIABLES 

102 JOPTIMsl 
INC3=1 

1 F MPHF *DK 2FK( TEMP (1,1)) 

TEMPCFsl)K2FK (T£MP(9, 1) ) 

IRITE*A 

C002*0R(8) /RAD(9) 

C0U=.5*C0D2 

OR 8 sO, 

IF (PEREO.GT.O..AND.VREQ.GT.O,) GO TO 105 
I R I R N a 3 
WRITE (6,108) 

10 b FORMAT ( 1 H// 100 H VOLTAGE AND POWER MUST BE SPECIFIED IN PUMP MOOUL 
IE PARAMETRIC CALCULATIONS, CALCULATIONS SUPPRESSED) 

return 

10b CUREQspEREQ/VREU 
RLPC ( .1) sVREQ/CUREO 
I Z9*2 

Z9SAV*DR(5) 

P»R01NC*FL0AT (NPUMP) 

R9MAX*Z9SAV+0,5*P 
0R3SAV*DR (3) 

ORbSAVsDR ( 6 ) 

DR3MAXs0R3SAV*P 

DRbMAXsDRbSAV+P 

103 WRITE (6,100) (CMT(J),J«1,25) 

100 FORMAT (lHl,ax,2A4, IX, 18AA,2H ( , F 2 . 0 , 2 ( 1 H/ , F 2 , 0 ) , 2 H) ,?A«) 

WRITE (6,101) TEMPHF,TEMKF,TEMPCF,TEMKF»ZID,0RTE,0R(1) ,pcmult,cod 

101 FORMAT (26H AVEHAGF HOT CLAD TEMP. *,FB.1»6H DEG. , A 1 , 3QX , 2 bHA V£R 
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1 AGE COLO CLAD TEMP, p,F 8.1*6H OEG, , Al / 1 OX , 1 fi>H INNER DIAMETER «, 

2 F8.4,4H IN, , 3bX , 22HT /E RADIAL THICKNESS 8,F8.4,4H IN./26H INNE 
3R CLAD THICKNESS s,F8,4,4H IN, , 8X , lbHNO. OF COUPLES «,F3.0,4X, 

4 27HCLA0 THICKNESS/O.D, RATIO «,F8,5) 

J»1 

IF (WTcON.GT.O.) Ja 2 

WRITE (6,104) CRMD(1,IENDK),CRmD(1,IENDK),CRMD(1,J) 

104 FORMAT (95H0 DR(3) QR(b) 0,0. ZLN ZLP RGEN 

1 VOC VLOAO CURRENT POWER 0 A4 , 5X , 3HET A A2 , 3X , A4 , 3HWT , 

2/119H (IN,) (IN,) (IN.) (IN.) (IN.) (M-OHMS) (VOLT 

3S) (VOLTS) (AMPS) (MATTS) (WATTS) (PCT.) (LBS.)) 

201 JOPTIMsJOPTIM+1 
IRTRNsl 

IF (ABS(CUR£0«CURU))/CUREQ.LE.T0LTEM,0R.J0PTIM,GT.10) GO TO 210 

RCPC*(VQC(1)«VREQ/PCMULT) /CUREO-RPPC 

RCONOsRPCd )«*RPPC-RPN 

P*RCPC*RCPC*4,*RPN*RCUND 

IF (P.LT.O.) GO TO 209 

ZLNs.5*ZLN*(RCPC-SQRT(P))/RC0ND 

RETURN 

209 ZLN«ZLN*SORT(RPN/RcONO) 

JOPT IMs 10 

RETURN 

210 OHMaPCMULT*RPC(i) 

WTsWATe(P) 

0CV«PCMULT*V0C(1) 

VLD=PCMULT*VPC(1) 

P*VLD*CUR(1) 

Q*PCMUlT*QT (1)+ZKEND*DTM0D(1) 

E * P / Q 

IF (INC3.EQ.0) GO TO 212 

WRITE (6,211) 0R(3) ,DR(6),Z0D,ZLN,ZLP,0HM,0CV, VLD,CUR(1) ,P,Q,E,WT 

211 FORMAT (2H0 , 3F8 , 4 , 2F9 . 4, 3PF9. 4, 0P2F I 0. 4 , F9, 2, F9, 3, F9. 1 , 2PF 1 0,4 , 

1 0PF9.3) 

INC3*0 
GO TO 214 

212 WRITE (6,213) DR (6) , ZOD , ZLN, ZLP # OHM, OC V » VLD , CUR ( l) , P , Q f E , WT 

213 FORMAT (10X,2F8.4,2F9.4,3PF9.4,0P2F10,4,F9.2,F9.3,F9,i,2PF10.4, 

1 0PF9.3) 

C - INCREMENT OUTER CONDUCTOR RADIAL THICKNESS 

214 J OP TIM? 2 
IRTRN*2 

DR(6)«DR(6)4R0INC 

IF (DR(b) .LE.DRbMAX) RETURN 
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DW (b) snRfoSAV 

C- INCREMENT INNER CONDUCTOR RADIAL THICKNESS 
OR(3)s D R( 3)+ROINC 
IF (DR(3) .LE.DR3MAX) RETURN 
JQPTIMsl 
DR (3) s L) R 3 S A V 

C- INCREMENT T/E WASHER RADIAL THICKNESS 
Z9SAVsz9SAV+0.25*R0INC 
DRC9J sDRta)*0.25*R0lNC 
IF (Z9SAV.LE.R9MAX) RETURN 
IRTRN*3 
RETURN 
END 


FUNCTION WATE(WTPC) 

CUMMQN /MDCPOP/RAD(9),0R(fl)#TEMP(9,30O) #V0C(300) , RPC (300) , PE (300) , 

1 QT (300) »N0PTIM» VPCI300) , RLPC ( 300 ) , DTMOD (300) , Cl>R ( 300 ) , I.NR AD ( 9 ) , 

2 ZLP,Zl N, ZLI#ZLNLP,ZLPNI, ZLTg , n&T 1 , NGT2 r I TPER T , PCMULT , NODUMP , I VD , 

3 QGEN,QGENL,QTE»OOUT,ZID# Z0D#ZkEND, ZK9,ZKR,RPPC,RPN, toltem,trej» 

4 VREQ,pEREG,ZLREQ,WTCON,ZLPN, I nUMZ,RN,DGRF,DGTRF,DGLRF,RAPK, JDUMP 

5 DP(30),NC 

COMMON /M0DCPL/CURM0DC15) , IT(300) * DT TE , CONRN, CONRP , ZL I P , ZNC , 

1 ALR64,MAXTEM,0GAM,CGAM, INDI f RlNT4,RINTb,RC0Na,RC0N6 f C1FEH,C1FEC, 

2 CONST 1 r CONST 2, CONST 3, C0NST9, CONS T5,NR50PT,SEGTE, ICRMDL#R ADC 
COMMON /F0WT/DEN(8) #DENl r DEN3N r DEN2N,DEN3PrDEN2P,DCRH,DCRC 
ZLMsZLI ♦ *001 

PIL=3,i4i59*ZLTE 

DENfiJe (ZLPNI*(JCRH + ZLM*D£NI J/ZlTE 
DENC4) = CZLN*DEn3n + ZLP*DEn3P + 2.*ZLM*DENI J/ZLTE 
DENIS) s(ZLN*DEN2N+ZLP*DEN2P+2.*ZLM*dENU/ZLTE 
DEN(b) s (ZLPNI*DCRC+ZLM*DENI)/Z lTE 
CALCULATE WEIGHT of one COUPLE 
W TPCsQ . 

DO 10 I * 1 » 0 

10 WTPCaWTPC+PIL*DR(I ) *(RAD(I)+RAD(I*1) 

CALCULATE TOTAL MODULE WEIGHT 

WATE*WTCON<r (RADC9) ♦RAD(l) ) *(RAD(R)-R 
return 
END 

FUNCTION DK2FKCT) 

COMMON /MOOOPT/lPDWT, DR4 , DR I NcR » IR I 
1 RIINC, R0lNC,RiMAX r R9MAX r DIMIN, QREU, 

2,NPERT, ZNPERT#CKFP,CKFT,CTE0C»CTEIC» 

CONVERTS KELVIN TO FAHRENHEIT DEGREES (IF 
DK2FKsT*CKFT-CKFP 

return 

END 


) *DEn ( I ) 

ADC1))+PCMULT*WTPC 


TE,DRTE» JOPTIM r IZ9, IENDK, IPUN, 

thing, tcinc, tcmax, thmax 

DELTL»RLPERT,RLMAX,ZNMIN,Z9SAV 

NZ(O)sl) 
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FUNCTION SI (XTBL# YT8L# XX,NN, IND# INDLE# INDUE) 

DIMENSION XT8L(2) , YTBL (2) 

XsXX 
N«NN 

IND a INDICATES T YPE -OF -E X TR aPOL A T I ON THAT WAS USED (IF ANY) 
(INDaO INDICATES NO EXTRAPOLATION WAS NEEDED ON X) 

(IND*1 INDICATES LOWER EXTRAPOLATION WAS NEEDED ON X) 
(IN0s2 INDICATES UPPER EXTRAPOLATION WAS NEEDED ON X) 

INDLE « INDICATES T YPE-OF-LOWER-EXTR APOLAT ION TO BE USED ON X 
(INDLEal INDICATES LOWER EXTRAP. ON X IS TO BE LINEAR) 

(INDLEs2 INDICATES LOWER EXTRAP. ON X IS TO BE PARABOLIC) 

( I NDLE*3 INDICATES LOWER-LINEAR EXTRAP. ON X AND ERROR PRI 
CINDLE«4 INDICATES LOWER-PARABOLIC EXTR. ON X AND ERROR PR 
ClNDLEaS INDICATES LOWER EXTRAP. ON X IS TO BE FIRST TABlF 
INDUE a INDICATES TYPE-OF-UPPER-EXTR APOLAT ION TO BE USED ON X 

( indue* i indicates upper extrap, on x is to be linear) 

(IN0UEa2 INDICATES UPPER EXTRAP. ON X IS TO BE PARABOLIC) 

CINDUEc3 INDICATES UPPER-LINEAR EXTRAP. ON X AND ERROR PRI 
(INDUEb 4 INDICATES UPPER-PARABOLIC EXTR, ON X AND ERROR PR 
(iNDUEaS INDICATES UPPgR EXTRAP, ON X IS TO BE LAST TABLE 
XTBL a NAME OF INDEPENDENT VARIABLE TABLES 
YTBL * NAME OF DEPENDENT VARIABLE TABLES 
N s NUMBER-OF-POINTS IN EACH TABLE 

X a PARTICULAR VALUE OF INDEPENDENT VARIABLE 

SIsYTBl(I) 

INDaO 

IFCNN.LE.l) return 
C CHECK TO SEE IF LOWER OUT-OF-RaNGE EXTRAPOLATION WILL BE NEEDED 

IF (X-XTBL(l)) 120# 130# 150 

C LOWER OUT-OF-RANGE EXTRAPOLATION WAS FOUND NECESSARY (SET INDal) 

120 INDal 

IF (INDUE .EQ. S) RETURN 
130 11*2 

GO TO 254 

C CHECK TO SEE IF UPPER OUT-OF-RaNGE EXTRAPOLATION WILL BE NEEDED 
150 IF(X,LT.XTBL(N) ) GO TO 210 

C UPPER OUT-OF-RANGE EXTRAPOLATION WAS FOUND NECESSARY (SET IND*2) 

IND«2 
II*N 

GO TO (254# IS0#254, 180, 131) # INDUE 
131 SlaYTBL(N) 

RETURN 
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GO TO 254 
C 

C (X IS IN-RANGE hake A POINT SEARCH On X TO OBTAIN IT) 

210 NMl*N-i 

00 220 IKs2,NMl 

1 I s I K 

IF (X TBlC IK )-X) 220,254,254 
220 CONTINUE 

c 

254 XlsXTBLUlM) 
x2=XTBLtII) 

YUYTBl(II-I) 

Y2*YTBlC1I) 

C 

C CHECK IF CUPPER OR LOWER) EXTRAPOLATION WAS FOUND TO Bp NECESSARY 

IF(IND-1)H59,257,258 
C 

C LOWER EXTRAPOLATION IS NEEDED - CHECK I NQLE FOR TYPE To BE USED 

257 GO TO C270, 259, 370,359), INDLE * 

C 

C UPPER EXTRAPOLATION IS NEEOED - CHECK INDUE FOR TYPF Tn BE USED 

258 GO TO <270,259,370,359) , INDUE 
C 

C ERROR PRINTOUT 

359 CALL ERROR ( 33H TABLE EXTRAPOLATED PARABOl ICALL Y) 

C 

259 IF (NN.LE.2) GO TO 270 
X3=XT8l (IIM) 

Y3=YTBL(II+1) 

C 

C' PARABOLIC INTERPOLATION OR EXTRAPOLATION 
2bO SI«Ylf(l,*(X2-X)/(X3-Xl) ) * { Y 2-Y 1 ) * C X-X 1 ) / ( X2-X 1) - C X2-X ) / ( X3-X 1 ) * ( 
lX-Xl)/fX3-X2)*(Y3-Y2) 

GO TO 300 
C 

C ERROR PRINTOUT 

3/0 CALL ERROR (28H TABLE EXTRAPOLATED LINEA«LY) 

C 

C LINEAR EXTRAPOLATION 

270 SI«Y14(Y2-Y1)*(X-X1) /(X2-X1) 

300 RETURN 

end 
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